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Abstract 

The  Hansen  Solubility  Parameters  (HSP)  for  several  Polyhedral  Oligomeric  SilSesquioxane  (POSS) 
compounds  were  successfully  determined,  demonstrating  the  applicability  of  the  HSP  approach  for 
selected  types  of  organic-inorganic  compounds.  As  commonly  practiced  with  organic  polymers,  a  set  of 
simple  “pass/fail”  tests  for  complete  solubility  at  a  fixed  concentration  (100  mg/mL)  was  conducted  for 
an  array  of  five  octameric  POSS  compounds,  octa(phenethyl),  octa(styrenyl),  octa(/.sobutyl), 
octakis(hexafluoroAobutyl),  and  (l-naphthyl)heptaphenyl),  and  45  test  solvents.  Group  contributions 
for  the  octameric  POSS  cage  were  determined  using  three  different  approaches,  which  produced  similar 
results.  The  best  cage  contribution  estimate  for  the  dispersive,  polar,  and  hydrogen-bonding 
components  8d,  Sp,,  and  5h  of  the  total  solubility  parameter  was  determined  to  be  8d  =  22  (J/cc)  ,  8p  = 
19  (J/cc)  ',  and  8h  =  15  (J/cc)  ,  with  an  estimated  uncertainty  of  approx.  5  (J/cc)  .  The  utility  of  the 
HSP  approach  was  demonstrated  by  successfully  identifying  mixtures  of  poor  solvents  that  provided 
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significantly  enhanced  solubility  for  octa(z'sobutyl)  POSS,  and  by  successfully  estimating  the  HSP  of 
octakis(trifluoropropyl)  POSS  from  group  contributions  derived  solely  from  aromatic  POSS 
compounds. 

1.  Introduction 

Over  the  past  few  decades,  Hansen  Solubility  Parameters  (HSP)  have  achieved  a  high  level  of 
prominence  in  the  coatings  industry  as  a  practical  tool  for  estimating  a  variety  of  thermodynamic  and 
transport  properties  in  polymer  systems.1’4  An  advantage  of  the  HSP  approach  over  the  simpler  one- 
component  Hildebrand  solubility  parameter  theory  is  that  the  HSP  approach  enables  the  identification  of 
mixtures  of  non-solvents  that,  when  combined  in  the  proper  ratio,  become  good  solvents  for  difficult  to 
dissolve  polymers.5  The  adoption  of  the  HSP  approach  to  the  more  general  problem  of  finding  sets  of 
polymers  and/or  small-molecule  fluids  that  display  desired  thermodynamic  interaction  characteristics, 
whether  favorable  or  unfavorable,  has  saved  countless  hours  of  trial  and  error  in  gas  separation 
membranes,6  drug  delivery,7’9  and  nanocomposite10’15  applications.  It  is,  therefore,  virtually  certain  that 
successful  efforts  aimed  at  expanding  the  range  of  applicability  of  Hansen  Solubility  Parameters  will 
generate  numerous  very  substantial  positive  impacts  across  a  wide  variety  of  industries. 

Polyhedral  Oligomeric  SilSesquioxanes  (POSS)16  are  a  family  of  inorganic/organic  core/shell 
nanostructures  possessing  molecular  weight  values  on  the  order  of  1000  g/mol,  in  which  the  core 
consists  of  a  silsesquioxane  cage  with  the  formula  (SiOi.5)n,  where  n  is  generally  between  8  and  14.  The 
shell  consists  of  n  organic  functional  groups  originating  at  each  Si  atom  on  a  cage  vertex.  Cages  are 
generally  categorized  according  to  the  bonding  type  and  number  of  silicon  atoms  that  comprise  each 
vertex.  The  functionality  may  vary  from  one  organic  group  to  another,  although  the  most  commonly 
used  compounds  feature  octameric  cages  with  identical,  simple  functional  groups. 

Since  their  commercialization  in  the  late  1990s,  POSS  compounds  have  found  a  large  variety  of  uses, 
particularly  as  modifiers  for  polymers  and  inorganic  fillers  in  the  medical,  aerospace,  and  electronics 
industries,19  where  they  serve  as  processing  aids  for  high-temperature  thermoplastics,  dispersion  and 
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compatibilization  agents,  and  thennal  and  electrical  insulation  enhancers.  The  inorganic  core  is  both 
mechanically  robust,  resistant  to  oxidation,  and  thermally  stable,  and  the  ability  to  synthesize  a  variety 
of  peripheral  groups  makes  the  tailoring  of  properties  fairly  straightforward.  Despite  the  many 

applications  for  POSS  compounds  that  involve  thermodynamic  compatibility  and  mixing,  there  have 

20  2 1 

been  few  systematic  comparative  attempts  to  quantify  POSS  solubility  in  common  organic  solvents.  ’ 

More  significant  effort,  however,  has  been  devoted  to  exploring  the  solubility  of  inert  POSS 
compounds  in  polymers.16, 19  There  have  been  numerous  reports  of  good  compatibility  between  certain 
POSS  compounds  and  selected  polymers,  e.g.,  octakis(2-phenylethyl),  referred  to  as  octa(phenethyl), 
POSS  and  either  polystyrene  or  polyvinyl  chloride.  However,  in  many  cases,  the  reported  solubility 
of  POSS  compounds  in  polymers  of  interest  is  quite  limited.  Therefore,  there  is  a  clear  need  for  a 
rational  approach  to  determine  the  solubility  and  thermodynamic  interactions  of  POSS  compounds  with 
both  solvents  and  polymers. 

Despite  the  obvious  need,  there  has  been  only  a  small  amount  of  published  work  to  date  on  the 
solubility  parameters  for  POSS  compounds.  Recent  work  reported  by  the  Morgan  group24'26  has 
focused  on  determination  of  the  Hildebrand  solubility  parameters  of  POSS  compounds  with  partial 
success  in  predicting  polymer  compatibility.  Other  recent  work  by  Lim  et  al.  ’  has  shown  that 
calculation  of  the  HSP  for  the  peripheral  arms  of  POSS  compounds  (ignoring  the  cage  entirely)  provides 
limited  comparative  insight  into  polymer  compatibility  as  well.  In  most  cases,  the  Hansen  Solubility 
Parameter  approach  has  been  limited  to  purely  organic  compounds,  with  only  a  few  investigations 
perfonned  for  inorganic  groups  ’  (mainly  in  the  context  of  nanoparticles  or  surfaces  ’  ),  or  for 

organic/inorganic  polymers,  such  as  polysiloxanes  '  and  polyphosphazines.  The  size  and  even 
shape  of  POSS  molecules,  however,  would  not  be  expected  to  preclude  the  determination  of  HSP 
because  the  HSP  approach  has  been  successfully  applied  to  compounds  such  as  fullerenes,12,  40 
asphaltenes,41'44  and  bitumen45,  46  that  are  intermediates  between  small  molecules  and  polymer  chains. 
Though  for  nanoparticles  in  general,  one  would  expect  the  solubility  parameters  to  reflect  the  surface, 
rather  than  bulk  composition,  for  POSS  compounds,  the  central  cage  is  at  least  partially  accessible  to 
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surrounding  molecules,  and  the  overall  issue  of  solvent  accessibility  is  similar  to  those  encountered  for 
many  types  of  polymers.  In  short,  given  the  many  benefits  that  will  accrue  from  successful  application 
of  HSP  to  POSS  compounds,  the  subject  is  worth  considerable  further  investigation. 

Herein  we  report  the  results  of  our  efforts  to  determine  Hansen  Solubility  Parameters  for  selected 
POSS  compounds.  We  found  that,  for  an  initial  group  of  five  POSS  compounds,  simple  solubility 
testing  of  the  type  carried  out  for  polymers  was  adequate  for  the  determination  of  the  HSP  values.  We 
further  found  that  the  consistency  of  the  results  was  of  the  same  quality  encountered  for  most  polymer 
systems,  and  that  there  appeared  to  be  no  fundamental  disadvantage  associated  with  POSS  compounds 
as  a  group.  Additional  analysis  of  these  HSP  values  allowed  us  to  investigate  the  possibility  of 
assigning  a  group  contribution  to  the  octameric  cage,  where  we  determined  that,  although  a  reasonably 
consistent  set  of  parameters  could  be  obtained  for  many  compounds,  the  issue  of  cage  accessibility  may 
complicate  matters  more  than  in  a  typical  polymer.  Our  results  represent  an  important  extension  of  HSP 
in  the  realm  of  inorganic/organic  materials,  a  significant  but  logical  step  forward  that  is  certain  to  bring 
substantial  benefits  to  the  ever  growing  array  of  applications  in  which  POSS  compounds  are  featured 
prominently. 

2.  Experimental 

2.1  Materials.  Octa(phenethyl)  POSS,  octa(styrenyl)  POSS,  octakis(hexafluoroAobutyl)  POSS,  and 
(l-naphthyl)heptaphenyl  POSS  were  synthesized  at  AFRL.  Octa(Aobutyl)  POSS  was  obtained  from 
Hybrid  Plastics  (Hattiesburg,  MS).  The  chemical  structure  of  these  compounds  is  given  in  Figure  1.  A 
large  amount  of  data  is  available  for  the  octa(Aobutyl)  compound,16  with  limited  data  available  for  the 
octa(phencthyl),'  ’  ’  '  octa(styrenyl),  ’  octakis(hexafluoroAobutyl),  and,  more  recently,  the  (1- 
naphthyl)heptaphenyl  compounds,  respectively.  The  set  of  45  solvents  used  in  the  study  were 
obtained  from  a  variety  of  commercial  sources.  Most  were  of  technical  grade,  and  all  were  used  in  as- 
received  condition.  All  are  listed  in  Supporting  Information. 
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la:  R,  R’ = -CH2CH2-C6H5 
1b:  R,R'  =  -CH=CH-C6H5 

Id:  R,R' = -CH2(CH)(CF3)2 
1e:  R  =  -C6H5;  R'  =  1 -naphthyl 

Figure  1 :  Chemical  structure  of  POSS  compounds  for  which  HSP  were  determined. 

2.2  Test  Procedures.  For  each  of  the  tests  reported  herein,  50  mg  of  the  POSS  compound  was  placed 
in  a  clean  5  mL  glass  vial  along  with  0.5  mL  of  the  selected  solvent.  The  tightly  closed  vials  were  then 
stirred  for  two  minutes  and  checked  for  dissolution.  When  complete  dissolution  did  not  occur,  the  vials 
were  allowed  to  stand  for  up  to  one  hour.  The  tests  were  then  rated  as  follows:  “Pass”  indicated 
complete  dissolution  within  two  minutes;  “Pass  — “  indicated  complete  dissolution  within  one  hour  but 
not  within  two  minutes  (in  these  cases,  most  of  the  solute  dissolved  initially  and  extra  time  was  required 
to  achieve  complete  dissolution,  an  indication  of  near  saturation  given  the  low  molecular  weight  of  the 
POSS  compounds);  “Fail  +”  indicated  that  significant  but  incomplete  dissolution  occurred;  “Fail” 
indicated  that  little  or  no  dissolution  was  observed.  Note  that,  in  general,  “Pass  — “  was  treated  as  “Pass” 
and  “Fail  +”  as  “Fail”. 

Analysis  Procedures.  The  primary  goal  of  the  analysis  for  a  given  POSS  compound  was  to  provide  a 
best  estimate  of  the  dispersive  (5D),  polar  (5P),  and  hydrogen  bonding  (8H)  components  of  the  Hansen 
Solubility  Parameter.  In  addition,  some  indication  of  the  uncertainty  in  the  HSP  values,  as  well  as  an 
estimate  of  the  “radius  of  interaction”  (Ro)  that  defines  the  expected  boundary  in  “HSP  space”  between 
“good”  and  “poor”  solvents,  was  also  desired.  More  information  on  these  concepts  can  be  found  in  the 
work  by  Hansen. 1 

The  basic  approach  was  to  maximize  the  goodness  of  fit  of  the  data  set  to  the  “sphere  of  solubility” 
described  by  Hansen,1  based  on  the  pass/fail  ratings  (here  denoted  as  “test;”  and  the  associated  HSP  (8D- 
s,i,  8P_s,i,  5H-s,i)  of  each  solvent  tested.  Namely,  for  each  instance  i,  a  characteristic  distance  Ra  was 
defined  such  that 


ic:  R,R'  =  -CH2(CH)(CH3)2 
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Ra,i—  4  *  (8d-P0SS  -  8o-S,i)"  +  (8p-P0SS  -  8p-s,i)“  +  (8h-P0SS  -  8||-s.,)2 


(1) 


with  the  criteria  being 

If  Ra,i<  Ro,  then  test;  =  “Pass.”  Otherwise,  testi  =  “Fail” 
in  which  8d-poss,  8p_poss,  8h-poss,  and  Ro  are  adjustable  parameters.  The  goodness  of  fit  was 
determined  by  means  of  a  x  test,  using  the  total  fraction  of  tests  grades  as  “Pass”  and  the  total  fraction 
of  solvents  for  which  Rajj  <  Ro  to  determine  the  distribution  of  expected  occurrences.  The  procedure 
described  above  was  the  computational  equivalent  of  attempting  to  identify  the  central  coordinates  and 
radius  of  a  sphere  that  included  all  solvents  with  a  “pass”  rating  but  excluded  all  solvents  with  a  “fail” 
rating,  in  a  Cartesian  space  with  axes  given  by  2  8d,  8p,  and  8h- 
To  locate  the  optimal  values  of  8d-poss,  8p_poss,  8h-poss,  and  Ro,  a  brute  force  method  based  on  an 
iteratively  refined  search  space  was  employed.  In  the  first  iteration,  Sd-poss  was  varied  from  15  to  20  in 
increments  of  0.5,  while  Sp.poss  and  8h-poss  were  varied  from  0  to  10  in  increments  of  1,  while  Ro  was 
varied  from  0  to  20  in  increments  of  0.25.  The  initial  range  represented  a  compromise  between 
comprehensiveness  and  the  need  to  minimize  computation  time,  and  was  based  on  the  range  of  the  HSP 
noted  for  solvents  with  a  rating  of  “pass”.  For  every  possible  combination  of  the  four  parameters  given 
the  ranges  and  increments  listed  above,  the  x  parameter  was  then  calculated  and  the  absolute  maximum 
recorded.  For  the  second  iteration,  the  search  grid  was  centered  on  the  point  having  maximum  x  (with 
the  lowest  value  of  all  other  parameters  used  to  break  ties)  found  in  the  previous  iteration,  with  the 
increments  reduced  to  0.4  for  Sd-poss,  and  0.8  for  Sp.poss  and  Sh-poss,  while  the  increment  was 
maintained  at  0.25  for  Ro.  (These  values  were  chosen  to  provide  conveniently  divisible  intervals  for 
subsequent  iterations  while  enabling  the  search  space  to  be  expanded  well  beyond  the  originally  chosen 
boundaries  if  necessary.)  The  process  was  repeated  two  additional  times,  centering  the  grid  as  described 
previously  and  reducing  the  increment  for  8d-poss,  8p_poss,  and  8h-poss  by  a  factor  of  4  each  time,  while 
maintaining  the  increment  for  Ro.  The  method  provided  a  final  search  resolution  of  0.025  for  8d-poss, 
0.05  for  Sp.poss  and  8h-poss,  and  0.25  for  Ro.  Note  that  the  finer  spacing  for  Sd-poss  was  used  to 
compensate  for  the  factor  of  4  in  the  Hansen  expression  (effectively,  the  parameter  2  8d-poss  was  used  as 
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a  dimension  along  with  Sp.poss,  an(l  §h-poss  in  a  uniform  Cartesian  search  grid  that  shrank  with  each 
iteration  until  the  optimal  center  co-ordinates  were  identified). 

Due  to  the  highly  non-linear  nature  of  the  optimization  process,  many  possible  locations  in  “HSP 
space”  bounded  by  an  irregular  and  even  non-contiguous  region  can  be  of  equal  or  nearly  equal 
goodness  of  fit,  thus  it  is  important  to  obtain  some  estimate  of  uncertainty.  To  address  this  issue,  we 
again  employed  a  brute  force  method,  calculating  the  %  value  for  all  combinations  of  8d-poss,  Sp-poss, 
5h-poss,  and  Ro  with  increments  of  0.05,  0.1,  0.1,  and  0.25,  respectively,  around  the  optimal  point 
identified  in  the  final  search  iteration  with  a  span  of  1 1  increments  for  8d-poss,  §p-poss,  and  8h-poss  and 
41  increments  for  Ro.  We  considered  %  values  that  differed  from  the  previously  determined  optimum  by 
less  than  one  to  be  representative  of  equivalent  goodness  of  fit.  This  practice,  originally  suggested  by  E. 
von  Meerwaal  of  the  University  of  Akron,  was  found  to  work  quite  well  in  this  situation.  The  values  of 
every  combination  of  5d-poss,  §p-poss,  and  8h-poss  for  which  at  least  one  value  of  Ro  enabled  the 
condition  for  equivalent  goodness  of  fit  to  be  met  were  logged.  If  any  of  these  coordinates  fell  on  the 
boundary  of  the  search  region,  the  increments  of  8d-poss,  §p-poss,  and  5h-poss  were  increased  in  linear 
steps  of  100%  (i.e.  they  were  first  doubled  from  their  original  values,  then  tripled  from  their  original 
values,  etc.),  and  the  entire  process  repeated,  until  no  co-ordinates  on  the  boundary  of  the  search  region 
met  the  criteria.  We  then  calculated  the  centroid  and  radii  of  gyration  (in  each  dimension)  of  the  set  of 
all  co-ordinates  that  met  the  goodness  of  fit  criteria,  using  the  centroid  for  the  finally  determined  values 
of  the  HSP  for  the  given  POSS  compound.  The  Ro  values  from  the  list  of  all  co-ordinate  combinations 
that  met  the  goodness  of  fit  criteria  were  averaged  and  reported  as  the  final  Ro  value  for  the  POSS 
compound.  The  radii  of  gyration  were  then  reported  as  the  “characteristic  uncertainty”  of  the 
measurement. 

3.  Results  and  Discussion 

3.1  General  Overview 

Figures  2-6  depict  the  outcomes  of  the  solvent  testing  for  the  octa(phenethyl)  (Figure  2), 
octa(styrenyl)  (Figure  3),  octa(Aobutyl)  (Figure  4),  octakis(hexafhroroAobutyl)  (Figure  5)  and  (1- 
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naphthyl)heptaphenyl  (Figure  6)  POSS  compounds.  In  each  of  these  Figures,  the  axes  represent  the 
HSP  with  the  scaling  suggested  by  Hansen,  with  5d  compressed  by  two.  The  fdled  symbols  represent 
the  HSP  of  solvents  that  rated  “pass”,  while  the  unfilled  symbols  reflect  the  solvents  that  rated  “fail”. 
The  spherical  object  represents  the  predicted  “pass/fail”  boundary,  based  on  the  calculated  HSP  for  the 
individual  POSS  compound  (the  center  co-ordinates  of  the  sphere)  and  the  calculated  value  of  the 
“radius  of  interaction”  Ro.  Numerical  values  of  the  calculated  HSP,  Ro,  characteristic  uncertainties  (see 
Experimental  section  for  details),  and  data  on  the  overall  solubility  have  been  provided  in  Table  1,  as 
well  as  the  number  of  failed  predictions  based  on  the  relative  energy  difference  (RED)  concept,1  in 
which  the  RED  is  defined  as  the  ratio  of  the  distance  in  “HSP  space”  between  the  co-ordinates  of  the 
POSS  compound  and  the  solvent  to  Ro,  with  “good”  solvents  expected  when  RED  <  1  and  “poor” 
solvents  expected  for  RED  >  1 . 


X 

<0 


Figure  2.  Results  of  solubility  testing  for  octa(phenethyl)  POSS.  Filled  symbols  indicate  solvents  with 
a  rating  of  “pass”  (complete  solubility  at  100  mg/mL),  while  open  symbols  indicate  solvents  with  a 
rating  of  “fail”.  The  “sphere  of  solubility”  that  best  describes  the  data  (based  on  maximizing  y')  is  also 
shown. 
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Figure  3.  Results  of  solubility  testing  for  octa(styrenyl)  POSS,  with  data  displayed  as  in  Figure  2. 


Figure  4.  Results  of  solubility  testing  for  octa(/.sobutyl)  POSS,  with  data  displayed  as  in  Figure  2. 
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Figure  5.  Results  of  solubility  testing  for  octakis(hexafluorozsobutyl)  POSS,  with  data  displayed  as  in 
Figure  2. 


i 
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Figure  6.  Results  of  solubility  testing  for  (l-naphthyl)heptaphenyl  POSS,  with  data  displayed  as  in 
Figure  2. 
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Table  1 


Estimated  HSP  and  Selected  Related  Characteristics  of  Octameric  POSS  Compounds  (Standard 


Analysis) 


Name  (Code) 

8d 

Centroid 

5d 

Rg 

5P 

Centroid 

Sp 

Rg 

SH 

Centroid 

5h 

Rg 

Ro 

#  Good  / 

#  Tested 

#  Anomalies 

Octa(phenethyl) 

(PE-1) 

19.9 

0.4 

7.4 

0.5 

6.3 

0.6 

9.8 

29/45 

2 

Octa(styrenyl) 

(ST-1) 

20.2 

0.7 

5.8 

0.8 

6.6 

0.7 

9.8 

21/45 

5 

Octa(z'sobutyl) 

(IB-1) 

18.0 

0.1 

2.1 

0.2 

2.7 

0.3 

4.3 

8/45 

6 

Octakis(hexafluoro- 
/.vobutyl)  (HF-1) 

15.3 

0.8 

9.3 

0.9 

11.0 

0.6 

7.3 

10/36 

2 

( 1  -naphthyl)hepta- 
phenyl  (NP-1) 

16.9 

0.1 

4.0 

0.2 

6.3 

0.2 

2.5 

3/45 

1 

Rg  =  radius  of  gyration,  a  measure  of  the  uncertainty 


Perhaps  the  most  important  general  result  that  can  be  appreciated  from  examination  of  Figures  2-6  is 
the  level  of  success  of  the  HSP  approach  for  the  POSS  compounds.  In  general,  the  best  fit  predictions 
result  in  errors  (counting  both  anomalous  good  and  anomalous  poor  solvents)  of  only  approx.  25%,  as 
often  as  they  successfully  identify  good  solvents,  whereas  assigning  otherwise  randomly  chosen  guesses 
based  only  on  the  proportion  of  solvents  rated  “pass”  for  a  particular  compound  results  in  about  2.5 
errors  for  every  good  solvent  successfully  identified.  In  other  words,  the  use  of  HSP  for  POSS 
compounds  results  in  about  a  90%  reduction  in  time  lost  to  trial  and  error. 

3.2  HSP  of  Individual  POSS  Compounds 

Of  the  POSS  compounds  examined,  octa(phenethyl)  POSS  (Figure  2)  proved  to  be  the  most  generally 
soluble,  achieving  a  rating  of  “pass”  in  29  out  of  the  45  solvents  tested.  All  but  two  of  the  passing 
solvents,  and  none  of  the  failing  solvents,  lay  inside  a  sphere  of  radius  9.8  centered  at  (6d,  5p,  5h)  = 
(19.9,  7.4,  6.3).  The  two  exceptions  were  cyclohexane  and  diethyl  ether,  which  were  only  slightly 
outside  the  sphere.  The  uncertainty  in  the  determined  HSP  values  could  be  traced  to  a  lack  of  non¬ 
solvents  with  a  suitably  high  8d  value,  which  results  in  less  than  optimal  constraint  on  the  location  of  the 
sphere  of  solubility.  Such  a  lack  of  optimal  constraint  is  a  fairly  common  issue  when  HSP 
determinations  are  performed. 
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Two  additional  attempts  to  generate  an  octa(phenethyl)  POSS  solubility  data  set  that  was  more  easily 
constrained  were  made.  First,  the  set  of  solubility  tests  was  repeated  using  a  concentration  of  200 
mg/mL,  rather  than  100  mg/mL,  with  the  intent  of  generating  fewer  passing  solvents  and  shrinking  the 
radius  of  interaction.  Although  the  number  of  passing  solvents  was  reduced  to  19,  this  approach  was 
unsuccessful  because  the  solvents  with  the  highest  5d  still  passed.  By  disproportionately  eliminating 
solvents  that  helped  to  constrain  the  sphere,  the  uncertainties  were  in  fact  made  worse.  (Data  for  this 
and  other  “Supplemental”  tests  is  provided  in  Table  2).  Secondly,  an  additional  two  solvents  and  ten 
solvent  mixtures  (listed  in  Supporting  Information)  were  chosen  based  on  their  location  in  “uncharted” 
regions  of  HSP  space  near  the  probable  boundary  of  the  sphere  of  solubility  and  tested  at  100  mg/mL. 
The  results  were  then  combined  with  those  from  the  original  45  solvents.  The  second  procedure  was 
modestly  effective  at  further  constraining  the  sphere  of  solubility.  The  results,  as  seen  in  Table  2,  were 
changed  little  in  terms  of  actual  HSP  co-ordinates.  However,  the  uncertainties  were  reduced. 


Table  2 

Estimated  HSP  and  Selected  Related  Characteristics  of  Octameric  POSS  Compounds  (Supplemental 
Analysis) 


Analysis  (Code) 

§D 

Centroid 

§D 

Rg 

Sp 

Centroid 

sP 

Rg 

SH 

Centroid 

SH 

Rg 

Ro 

#  Good  / 

#  Tested 

# 

Anomalies 

Octa(phenethyl)  @  200 
mg/mL  (PE-A1) 

21.6 

0.5 

6.2 

0.2 

9.9 

0.4 

11.8 

19/45 

2 

Octa(phenethyl)  with 
additional  solvents  and 
mixtures  (PE-A2) 

19.7 

0.3 

8.0 

0.2 

5.6 

0.4 

9.6 

35/57 

3 

( 1  -naphthyl)hepta- 
phenyl  @  50  mg/mL 
(NP-A1) 

16.9 

0.1 

4.0 

0.2 

6.3 

0.2 

2.5 

3/45 

1 

( 1  -naphthyl)hepta- 
phenyl  with  “fail  +” 
considered  equal  to 
“pass”  (NP-A2) 

18.1 

0.2 

8.6 

0.3 

6.7 

0.1 

4.3 

7/45 

2 

The  solubility  of  octa(styrenyl)  POSS,  a  more  rigid  analog  of  octa(phenethyl)  POSS,  was  predicted 
only  modestly  less  accurately  by  the  “spherical  rule”,  with  five  exceptions  (all  non-solvents  within  the 
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sphere)  out  of  45  solvents,  and  with  three  of  the  exceptions  exhibiting  partial  solubility.  Though 
octa(styrenyl)  POSS  was  somewhat  less  soluble  overall  (21  passing  solvents)  compared  to 
octa(phenethyl)  POSS,  the  Ro  values  for  octa(styrenyl)  and  octa(phenethyl)  POSS  were  equivalent. 
The  octa(styrenyl)  POSS  had  a  slightly  higher  5d  value  than  octa(phenethyl)  POSS,  meaning  that  its 
calculated  HSP  also  were  made  more  uncertain  by  lack  of  constraints. 

Though  known  for  its  good  polymer  compatibility  characteristics,16  octa(Aobutyl)  POSS  (Figure  4) 
exhibited  a  significantly  smaller  radius  of  interaction  than  either  octa(phenethyl)  or  octa(styrenyl) 
POSS.  Though  well-constrained,  the  solubility  envelope  appeared  less  well-defined,  with  six 
exceptions  (three  of  each  type)  found.  Given  that  only  eight  solvents  passed  (with  eight  also  found 
inside  the  sphere),  the  odds  that  a  solvent  having  sufficiently  similar  HSP  to  octa(Aobutyl)  POSS  was, 
in  fact,  a  good  solvent  were  just  over  60%.  On  the  other  hand,  good  solvents  for  octa(Aobutyl)  POSS 
were  found  at  RED  values  (the  distance  to  the  center  of  the  sphere  divided  by  the  radius  of  interaction) 
of  1.6.  Interestingly,  the  anomalous  poor  solvents  (benzene,  xylenes,  and  phenyl  acetylene)  were  all 
small  aromatic  hydrocarbons,  while  toluene  was  soluble  at  100  mg/mL).  The  anomalous  good  solvents 
were  n-hexane,  diethyl  ether,  and  tetrahydrofuran.  While  the  predictive  power  for  the  HSP  approach 
was  less  effective  for  octa(Aobutyl)  POSS  than  for  any  other  compound  in  this  study,  it  was  still  on  par 
with  test  results  from  many  organic  polymers. 

In  order  to  further  explore  similar  POSS  compounds,  testing  was  undertaken  with  36  solvent 
candidates  on  octakis(hcxafluoro/.sobutyl)  POSS  (Figure  5),  a  member  of  the  class  of  compounds  known 
as  “fluoro  POSS”  that  have  played  a  key  role  in  realizing  superoleophobicity  on  rough  surfaces.55, 56  57 
The  octakis(hexafluoroAobutyl)  POSS  was  in  general  more  soluble  than  octa(Aobutyl)  POSS,  having  a 
larger  radius  of  interaction.  There  were  also  just  two  exceptions  to  the  “spherical  rule”,  /.sopropanol  (an 
anomalous  poor  solvent  lying  just  inside  the  boundary  of  the  sphere,  note  that  1 -propanol  and  1 -butanol 
were  good  solvents  and  within  the  sphere),  and  n-heptane  as  an  anomalous  good  solvent. 

Finally,  a  highly  aromatic  POSS  compound  was  also  examined,  (l-naphthyl)heptaphenyl  POSS 
(Figure  6).  Although  most  aromatic  POSS  compounds  are  very  poorly  soluble,  the  asymmetric 
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2 1 

substitution  pattern  appears  to  enable  at  least  limited  solubility.  A  total  of  three  good  solvents  enabled 
a  rudimentary  estimate  of  the  HSP  for  the  compound,  with  a  small  radius  of  interaction  and  just  one 
exception  to  the  “spherical  rule”,  albeit  a  crucial  one:  1 -methyl-2 -pyrrolidone  proved  to  be  a  good 
solvent  despite  having  an  RED  value  of  3.4.  With  only  three  good  solvents,  it  is  difficult  to  justify 
ignoring  this  anomaly  on  the  basis  of  only  two  other  data  points. 

In  an  effort  to  clarify  the  HSP  for  this  compound,  attempts  were  made  to  identify  more  “good” 
solvents  to  improve  the  robustness  of  the  analysis.  First,  the  non-solvents  were  all  re-tested  at  a  lower 
concentration  of  50  mg/mL  to  try  and  broaden  the  criteria  for  inclusion  of  “good”  solvents. 
Unfortunately,  in  every  case  for  which  the  aromatic  POSS  compound  failed  to  dissolve  at  100  mg/mL,  it 
also  failed  to  dissolve  completely  at  50  mg/mL.  However,  in  addition  to  the  three  “good”  solvents 
(chloroform,  tetrahydrofuran,  and  1 -methyl-2 -pyrrolidone),  there  were  four  solvents  in  these  tests  that 
received  the  “fail  +”  rating,  meaning  that  partial  solubility  was  noted.  Although  less  quantitative  than 
the  standard  approach,  adding  these  four  solvents  to  the  list  of  “good”  solvents  allowed  for  a  more 
robust  calculation  of  the  HSP,  with  the  results  shown  in  Table  2.  In  this  case,  no  poor  solvents  lay 
within  the  sphere  of  solubility,  and  only  two  good  solvents  (chloroform  and  diphenyl  ether)  were  found 
outside  of  the  sphere,  in  both  cases  by  a  short  distance. 

3.3  Estimation  of  Group  Contribution  Parameters  for  Octameric  POSS  Cages 

Having  identified  HSP  for  five  POSS  compounds,  it  is  reasonable  to  consider  whether  or  not  a 
systematic  pattern  exists  with  respect  to  the  HSP  of  the  peripheral  groups.  If  a  systematic  pattern  does 
exist,  then  it  should  be  possible  to  identify  a  “group  contribution”  for  the  octameric  POSS  cage  at  the 
core  of  the  five  compounds  studied.  With  such  a  “group  contribution”  in  hand,  it  would  then  be 
possible  to  calculate  estimated  HSP  for  any  octameric  POSS  compound.  (Other  methods  based  on  more 
sophisticated  computational  techniques58'64  might  also  be  useful  in  estimating  the  HSP  for  POSS 
compounds.  However,  a  group  contribution  approach  would  be  highly  valuable  due  to  its  simplicity.)  It 
should  be  noted,  however,  that  with  HSP  data  available  for  only  a  very  few  compounds,  any  such  group 
contribution  estimates  will,  of  necessity,  be  of  a  rough,  preliminary  nature. 
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In  order  to  determine  whether  a  meaningful  group  contribution  exists,  one  may  formally  divide  the 
cohesive  energy  of  the  POSS  compounds  into  two  parts,  one  arising  from  the  cage  and  one  arising  from 
the  peripheral  arms,  or 


Ei,P0SS  Ei  cage  T  Yij=i  ^i,jth-c 


(2) 


in  which  E;  represents  the  ith  cohesive  energy  component  (dispersive,  polar,  or  hydrogen  bonding)  for  a 
system  comprising  the  central  cage  and  n  peripheral  organic  groups  (arms),  not  all  of  which  are 
necessarily  identical.  This  formal  division  is  the  simplest  method  that  allows  for  modeling  of  a  wide 
variety  of  POSS  compounds.  If  values  for  Egth-arm  can  be  calculated  through  existing  group  contribution 
methods,  then  the  total  component  cohesive  energy  densities  for  the  POSS  cage  may  be  used  to  generate 
an  independent  estimate  of  all  Ei!Cage  for  every  POSS  compound  measured.  If  a  volume  contribution  for 
the  cage  can  be  determined,  then  component  solubility  parameters  for  the  cage  may  also  be  estimated. 
By  examining  the  consistency  of  the  values  of  EijCage  (or  the  cage  solubility  parameters)  derived  from 
different  compounds,  a  sense  may  be  developed  as  to  whether  a  single  meaningful  parameter  exists,  as 
well  as  a  sense  of  what  limitations  are  inherent  to  the  approach. 

In  Table  3,  the  experimentally  observed  volumetric  properties  of  POSS  compounds  are  compared 
against  values  estimated  by  three  different  methods:  1)  the  group  contribution  method  originally 
developed  by  Beerbower  for  solubility  parameter  estimation  as  expanded  and  reported  by  Hansen,1  2)  a 
more  recent  group  contribution  method  developed  by  Constantinou,65  and  3)  a  “best  available  analog 
group”  method  described  in  detail  in  Appendix  A.  Although  this  last  method  relies  more  heavily  on 
intuition  and  is  not  readily  generalized,  it  is  more  closely  linked  to  relevant  experimental  data.  The 
experimental  molar  volume  for  POSS  compounds  (obtained  by  gas  pycnometry  of  precipitated  powder, 
or,  in  the  case  of  octakis(hexafluoroAobutyl)  POSS,  X-ray  diffraction  data)  was  divided  into  a  portion 
representing  the  cage  (134  cc/mol,  based  on  a  model  described  in  Appendix  B)  and  a  remainder  that  is 
assigned  to  the  arms,  with  the  results  listed  in  Table  3.  Note  that  the  results  of  this  procedure  are  also 
displayed  for  octa(methyl)  POSS  as  a  further  check  on  the  validity  of  the  method. 
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Table  3 


Calculated  and  Experimental  Volumetric  Characteristics  of  POSS  Compounds 


Arm  Type 

Best  Available 
Analog  Group 
(BAAG)  relative 
weight) 

BAAG 
Volume 
(all  arms) 

Beerbower 
GC  Volume 
(all  arms) 

Constantinou 
GC  Volume 
(all  arms) 

Estimated  Arm 
Volume 
(experimental  - 
134  cc/mol) 

Experimental 

Volume 

(source) 

None 

Model 

134 

n/a 

n/a 

n/a 

n/a 

-CH3  x  8 

n/a 

n/a 

252 

209 

226 

360  (L) 

-Phenethyl  x  8 

Ethylbenzene 

845 

829 

868 

852 

986  (P) 

-Styrenyl  x  8 

Styrene 

785 

787 

799 

827 

961  (P) 

-/.vo  butyl  x  8 

/.so  butane 

639 

625 

617 

616 

750  (P) 

-Hexafluoro- 
/'.vobutyl  x  8 

Hexafluoro- 

zsobutylene* 

924 

870 

n/a 

743 

877  (X) 

-( 1  -naphthyl) 

+  -Phenyl  x  7 

Naphthalene  (1)  / 
Benzene  (7) 

642 

619 

619 

681 

815  (P) 

Given  the  many  assumptions  and  uncertainties  involved,  the  agreement  between  the  experimental  and 
theoretical  molar  volume  values  for  the  arms  shown  in  Table  3  is  quite  good,  with  the  exception  of 
octakis(hexafluoroAobutyl)  POSS.  In  that  case,  almost  all  of  the  methods,  including  the  experimental 
data,  are  less  reliable,  so  the  relative  lack  of  agreement  is  not  surprising.  Moreover,  no  single  method  is 
clearly  better  or  worse  than  the  others.  For  calculation  of  the  “group  contribution”  HSP  of  the  cage, 
however,  only  the  estimate  based  on  experimental  data  is  used  for  the  volume  of  the  cage  and  the  arms. 

Table  4  shows  estimates  for  the  cohesive  energy  components  of  the  arms,  again  based  on  three 
theoretical  methods,  the  best  available  analog  group  method  described  in  Appendix  A  (for  which  a 
direct  calculation  of  the  energy  components  is  employed),  the  Beerbower  group  contribution  method  as 
reported  by  Hansen1  (again,  for  which  a  direct  calculation  of  the  energy  components  is  available),  and  a 
more  recent  method  described  by  Stefanis.  Unlike  the  other  two  methods,  the  Stefanis  group 
contribution  involves  a  calculation  of  group  HSP,  these  are  then  multiplied  by  the  experimental  arm 
volume  to  obtain  the  dispersive,  polar,  and  hydrogen  bonding  cohesive  energy  components,  ED,  EP,  and 
Eh,  respectively.  The  Stefanis  method  is  available  without  modification  for  all  groups  needed. 
However,  the  Beerbower  method  is  not  directly  applicable  to  the  1 -naphthyl  group  or  the  fluorocarbon 
groups  for  all  components.  To  make  use  of  the  Beerbower  method,  two  assumptions  are  employed. 
First,  the  cohesive  energy  for  the  naphthyl  group  is  calculated  by  scaling  the  values  for  phenyl  groups 
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by  the  number  of  carbons  (10  for  naphthyl  as  opposed  to  6  for  phenyl);  second,  the  -CF3  group 
contribution  for  the  polar  component  is  calculated  by  taking  the  listed  value  for  the  total  energy  and 
subtracting  the  values  of  the  other  two  components  (which  are  listed  as  zero),  even  though  the 
contribution  of  the  polar  component  is  listed  as  “?”.  As  described  in  Appendix  A,  the  best  available 
analog  group  method  can  only  generate  an  estimate  for  the  dispersive  component  of  the  cohesive  energy 
of  the  hexafluoro/sobutyl  group.  Although  a  negligible  contribution  from  hydrogen  bonding  to  the 
cohesive  energy  density  of  this  group  represents  a  reasonable  assumption,  we  have  herein  also  used  a 
negligible  contribution  from  the  polar  component  for  lack  of  a  more  suitable  alternative. 

Table  4 


Calculated  Ann  Group  Cohesive  Energies  (J/cc)  for  POSS  Compounds 


Arm  Type 

Best  Avai 

able  Analog  Group 

Beerbower  GC 

Stefanis  GC 

Ed 

EP 

Eh 

Ed 

EP 

Eh 

Ed 

EP 

Eh 

-Phenethyl  x  8 

312000 

350 

1930 

331000 

1670 

1670 

320000 

1850 

10 

-Styrenyl  x  8 

320000 

930 

15500 

311000 

2880 

13700 

290000 

2440 

100 

-/.vo  butyl  x  8 

154000 

0 

0 

142000 

0 

0 

159000 

2360 

210 

-Hexafluoro- 
ivobutyl  x  8 

178000 

?(0) 

0 

90400 

67000 

0 

227000 

10500 

2490 

-(1 -Naphthyl) 

+  -Phenyl  x  7 

212000 

120 

2650 

273000 

1810 

1810 

246000 

2550 

210 

From  Table  4,  it  can  be  seen  that  the  agreement  among  the  different  methods  is  best  for  the  phenethyl 
and  zsobutyl  arms,  and  also  is  reasonable  for  the  polar  and  hydrogen  bonding  components  of  the  1- 
naphthyl/phenyl  arms,  along  with  the  hydrogen  bonding  component  of  the  hex  a  11  uo  ro/.vob utyl  arms.  For 
the  styrenyl  arm,  all  three  methods  are  in  reasonable  agreement  for  the  polar  component,  while  the  best 
available  analog  group  and  Beerbower  methods  agree  reasonably  well  for  the  dispersive  and  hydrogen 
bonding  components,  with  the  Stefanis  method  providing  somewhat  different  estimates. 

To  obtain  estimated  group  contribution  HSP  values  for  the  POSS  cages,  the  total  components  of  the 
cohesive  energy  were  calculated  based  on  the  experimental  solubility  parameters  and  molar  volumes, 
subtracted  from  the  calculated  values  for  the  arms,  and  then  divided  by  the  estimated  cage  volume  (134 
cc/mol,  from  Appendix  B)  before  raising  to  the  0.5  power.  These  group  contribution  component 
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parameters  are  shown  in  Table  5  for  the  best  available  analog  group  method,  Table  6  for  the  Beerbower 
method,  and  Table  7  for  the  Stefanis  method.  In  Tables  5-7,  a  measure  of  the  sensitivity  of  the 
calculated  HSP  group  contributions  has  also  been  provided.  The  sensitivity  estimate  reflects  the  effects 
of  assuming  an  error  in  the  measurement  of  HSP  equal  to  the  characteristic  uncertainty  (explained  in  the 
Experimental  section)  and  a  simultaneous  10%  error  in  the  estimate  of  ann  cohesive  energy,  with  the 
combination  always  done  in  such  a  way  as  to  augment  rather  than  offset  error.  For  example,  an 
overestimate  of  experimental  HSP,  and  thus  total  cohesive  energy,  is  combined  with  a  simultaneous 
underestimate  of  the  cohesive  energy  of  the  arms  to  maximize  the  remaining  cohesive  energy  assigned 
to  the  cage.  Note  that  calculations  based  on  the  alternate  analyses  of  the  experimental  data  are  also 
included  in  Tables  5-7. 

Considering  the  range  of  previously  measured  HSP  values,  it  is  reassuring  to  note  that,  although  in  a 
few  instances,  a  negative  cohesive  energy  contribution  for  the  cage  was  found,  in  no  case  is  the  resultant 
HSP  value  larger  than  30,  meaning  that  the  estimated  contributions  are  at  least  reasonable  in  magnitude. 
Moreover,  particularly  for  the  polar  and  hydrogen  bonding  components,  the  sensitivity  is  usually  small. 
For  the  dispersive  component,  the  sensitivity  values  are  significantly  larger.  The  greater  sensitivity  to 
error  results  from  the  fact  that  the  dispersive  energies  of  the  system  as  a  whole  are  not  much  different 
than  the  calculated  dispersive  energies  of  the  arms.  Therefore,  a  subtraction  of  two  large  numbers  is 
involved  to  determine  the  dispersive  component  contribution  of  the  cage.  This  fact,  in  and  of  itself, 
suggests  that  the  dispersive  component  of  the  cage  should  not  be  drastically  different  than  that  of  the 
arms,  while  the  polar  and  hydrogen  bonding  components  in  fact  should  be  quite  large  compared  to  the 
arms  in  all  cases  (with  the  possible  exception  of  the  polar  component  of  the  h e x a II  uo ro /.sob uty  1  arms). 
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Table  5 


Calculated  Group  Contribution  from  POSS  Cage  Based  on  HSP  Analysis  of  POSS  Compounds 
(Beerbower  Group  Contributions  for  Arms) 


Analysis 

Code 

Calculated 
Cage  8d 

Sensitivity 

Calculated 
Cage  8P 

Sensitivity 

Calculated 
Cage  8h 

Sensitivity 

PE-1 

21.0 

7.4 

19.8 

1.4 

16.7 

1.7 

PE-A2 

19.6 

7.8 

21.4 

1.4 

14.8 

1.8 

ST-1 

24.7 

7.7 

14.8 

2.3 

14.5 

2.6 

IB-1 

27.4 

2.2 

5.0 

0.7 

6.4 

0.7 

HF-1 

29.3 

3.7 

<0 

n/a 

28.1 

1.5 

NP-1 

<0 

n/a 

9.7 

1.9 

15.3 

1.1 

NP-A2 

<0 

n/a 

20.9 

1.8 

16.1 

1.0 

Table  6 

Calculated  Group  Contribution  from  POSS  Cage  Based  on  HSP  Analysis  of  POSS  Compounds 
(Stefanis  Group  Contributions  for  Anns) 


Analysis 

Code 

Calculated 
Cage  5d 

Sensitivity 

Calculated 
Cage  8P 

Sensitivity 

Calculated 
Cage  8h 

Sensitivity 

PE-1 

22.9 

6.8 

19.7 

1.4 

17.1 

1.6 

PE-A2 

21.6 

7.1 

21.4 

1.4 

15.2 

1.6 

ST-1 

27.6 

6.9 

15.0 

2.3 

17.6 

1.9 

IB-1 

25.1 

2.6 

2.7 

1.4 

6.2 

0.7 

HF-1 

<0 

n/a 

22.1 

2.6 

27.8 

1.6 

NP-1 

<0 

n/a 

8.9 

0.6 

15.5 

0.5 

NP-A2 

12.4 

7.1 

20.8 

0.8 

16.5 

0.2 

Table  7 

Calculated  Group  Contributions  from  POSS  Cage  Based  on  HSP  Analysis  of  POSS  Compounds  (Best 
Available  Analog  Group  Contributions  for  Arms) 


Analysis 

Code 

Calculated 
Cage  SD 

Sensitivity 

Calculated 
Cage  8P 

Sensitivity 

Calculated 
Cage  8h 

Sensitivity 

PE-1 

24.2 

6.4 

20.0 

1.4 

16.7 

1.7 

PE-A2 

22.9 

6.7 

21.6 

1.4 

14.7 

1.7 

ST-1 

23.1 

8.2 

15.3 

2.2 

14.0 

2.7 

IB-1 

25.8 

2.5 

5.0 

0.7 

6.4 

0.7 

HF-1 

14.2 

8.2 

23.8 

2.3 

28.1 

1.5 

NP-1 

12.3 

5.9 

9.8 

0.5 

14.9 

0.6 

NP-A2 

20.2 

4.5 

21.1 

0.7 

15.9 

0.3 
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As  expected,  extracting  a  meaningful  cage  contribution  from  the  data  on  octa(hexafluoroAobutyl) 
POSS  is  the  most  problematic.  The  main  issue,  however,  is  not  the  experimental  data,  but  rather  the 
inability  to  determine  a  reasonable  estimate  of  the  cohesive  energy  components  for  the  arms, 
particularly  for  the  dispersive  and  polar  components.  For  the  hydrogen  bonding  component,  the 
calculated  values  agree  reasonably  well  with  one  another,  but  are  considerably  higher  than  those  derived 
from  measurements  of  the  other  POSS  compounds.  One  possible  explanation  would  be  a  strong 
electronic  interaction,  as  suggested  by  X-ray  data,56  between  the  cage  and  the  fluorine  atoms  on  the 
arms,  at  least  some  of  which  must  approach  the  cage  closely  due  to  the  particular  geometry  of  the  arms. 

The  dispersive  and  polar  cage  contributions  derived  from  the  standard  analysis  of  the  (1-naphthyl)- 
heptaphenyl  POSS  also  appear  problematic.  However,  the  alternate  analysis  provides  values  that  are 
similar  both  among  methods  and  to  those  found  from  analysis  of  other  POSS  compounds.  In  the  case  of 
the  dispersive  component,  there  is  disagreement  among  the  computational  methods  as  to  the 
contribution  from  the  arms,  with  the  values  suggested  by  the  Beerbower  method,  and  to  a  lesser  extent, 
the  Stefanis  method  appearing  too  large.  In  this  case,  both  methods  have  known  shortcomings.  The 
Beerbower  method  does  not  provide  an  explicit  estimate  for  the  naphthyl  group,  while  the  Stefanis 
method  involves  an  explicit  dependence  of  the  cohesive  energy  density  on  molecular  size,  which  is 
likely  to  result  in  erroneous  values  for  larger  molecules  such  as  naphthalene.  As  for  the  polar 
component,  the  discrepancies  appear  to  result  principally  from  the  difficulties  with  the  standard  analysis 
(too  few  good  solvents  by  the  standard  criteria). 

A  final  discrepancy  is  observed  in  the  case  of  octa(/.sobutyl)  POSS.  Although  all  three  computational 
methods  provide  similar  results,  the  polar  and  hydrogen  bonding  components  inferred  from 
measurements  of  the  octa(zsobutyl)  POSS  are  much  lower  than  derived  from  the  other  POSS  cages. 
Among  the  POSS  compounds  analyzed,  octa(Aobutyl)  POSS  conformed  the  least  to  the  “spherical 
rule.”  However,  the  derived  cage  contributions  appear  relatively  insensitive  to  the  characteristic 
uncertainties  in  the  experimental  measurement.  It  appears  as  though  the  apparent  strongly  polar  and 
hydrogen  bonding  character  of  the  cage  is  screened  by  the  /.sobutyl  arms.  Such  an  effect  could  result 
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from  the  more  oblate  shape  of  the  Aobutyl  arms  in  comparison  to  the  other  POSS  types.  In  fact,  it  does 
seem  to  be  the  case,  based  on  the  limited  data  set  generated,  that  POSS  compounds  with  more  prolate 
arms  (such  as  phenethyl,  styrenyl,  and  (l-napthyl)heptaphenyl)  provided  more  consistent  cage 
contributions.  Such  a  “prolate  arm”  configuration  should  involve  less  direct  interaction  between  the 
cage  and  the  arms,  and  thus  ought  to  be  more  amenable  to  analysis  by  additive  group  contributions. 

Despite  all  of  the  difficulties,  it  does  appear  that,  at  least  for  the  prolate  arm  types,  a  “consensus” 
value  for  the  cage  contributions  is  available.  As  with  all  group  contribution  methods,  there  is  a 
considerable  amount  of  uncertainty.  However,  in  many  cases,  group  contribution  methods  are  utilized 
to  provide  rough  estimates.  In  these  cases,  the  values  provided  by  averaging  the  octa(phenethyl) 
(alternate  case  -A2,  which  appears  to  be  slightly  more  precise),  octa(styrenyl),  and  (1- 
napthyl)heptaphenyl  alternative  -A2)  would  seem  to  be  provide  the  best  estimate.  For  the  best 
available  analog  group  method,  the  averages  are  8d  =  22  (J/cm)  ,  8P  =  19  (J/cm)  ”,  and  Sh  =  15  (J/cm) 

1  /9 

,  or,  in  terms  of  cohesive  energy  components,  ED  =  65000  J/mol,  EP  =  48000  J/mol,  and  EH  =  30000 
J/mol.  Although  a  characteristic  uncertainty  is  not  readily  assigned  to  these  values,  based  on  the  data  in 
Tables  5-7,  most  of  the  values  for  the  three  compounds  lie  within  about  5  (J/cm)  of  these  averages. 
As  with  any  group  contribution  method,  HSP  data  for  many  additional  POSS  compounds  will  be  needed 
in  order  to  transform  the  rough  estimate  provide  above  into  a  highly  reliable  predictive  tool. 

3.4  Tests  of  the  HSP  and  Group  Contribution  Approaches  for  POSS  Compounds 

In  order  to  truly  evaluate  the  usefulness  of  the  HSP  approach,  as  well  as  the  group  contribution 
approach  to  estimating  HSP,  for  POSS  compounds,  two  specific  tests  of  predictive  power  were 
perfonned.  For  HSP  in  general,  one  of  the  most  useful  predictive  features  is  the  ability  to  identify 
mixtures  of  non-solvents  that  can  serve  as  a  good  solvent.  Thus,  guided  by  the  calculated  HSP  values, 
we  sought  to  identify  a  mixture  of  two  poor  solvents  that  displayed  anomalously  good  solubility  for  a 
POSS  compound.  For  this  demonstration,  we  chose  octa(/.vobutyl)  POSS  because  its  location  on  the 
HSP  diagram  and  its  relatively  small  radius  of  interaction  provide  for  many  possible  non-solvent 
combinations,  and  because  it  showed  a  large  number  of  exceptions  to  the  “spherical  rule”,  hence  it 
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serves  as  a  fairly  stringent  test  of  the  theory.  1,4-Dioxane  and  n-heptane,  both  of  which  are  poor 
solvents  for  octafi.sobutyl)  POSS,  were  chosen  for  the  experiment.  In  Figure  7,  quantitative  solubility 
data  for  octa(/.s'obutyl)  POSS  in  mixtures  of  1,4-dioxane  and  n-heptane  are  presented  as  a  function  of 
composition.  These  two  non-solvents  were  the  first  fully  miscible  pair  we  tested  (n-hexane  and  1,1, 2, 2- 
tetrabromoethane,  as  well  as  1-decene  and  quinoline,  which  we  also  attempted  to  test,  were  not  fully 
miscible  with  one  another). 
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Figure  7.  Solubility  of  octa(z'sobutyl)  POSS  at  20  °C  as  a  function  of  1,4-dioxane  fraction  in  a  mixed 
solvent  system  containing  n-heptane  and  1,4-dioxane. 


As  is  clearly  evident  in  Figure  7,  a  significant  improvement  in  solubility  is  obtained  when  n-hexane 
and  1,4-dioxane  are  utilized  as  a  mixed  solvent  system  for  octa(z.sobutyl)  POSS.  Although,  strictly 
speaking,  the  most  soluble  mixtures  fall  somewhat  short  of  the  “pass”  threshold  of  100  mg/mL,  despite 
the  fact  that  they  lie  just  within  the  sphere  of  interaction  on  the  HSP  diagram,  they  are  more  than  twice 
as  soluble  as  n-hexane  alone  and  about  four  times  as  soluble  as  a  simple  rule  of  mixtures  would  predict. 
Moreover,  the  most  soluble  mixtures,  containing  25-30  vol%  1,4-dioxane,  are  reasonably  close  in 
composition  to  the  value  (about  45  vol%  1,4-dioxane)  that  corresponds  to  the  minimum  RED  for  the 
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mixed  solvent  system  based  on  the  measured  HSP  values  for  octa(Aobutyl)  POSS.  Thus,  at  a 
qualitative  level,  the  predictions  of  HSP  theory  were  confirmed  for  octa(Aobutyl)  POSS. 

In  order  to  test  the  group  contribution  approach,  we  again  chose  a  stringent  challenge  for  the  theory. 

/TO 

The  compound  octakis(trifluoropropyl)  POSS  has  been  previously  reported  in  studies  of  liquid 
repellent  surfaces,69  and  is  potentially  useful  as  an  additive  to  fluoropolymers.  Therefore,  knowledge  of 
its  HSP  would  have  many  practical  benefits.  Currently,  however,  the  pure  compound  (containing  only 
octameric  cages)  is  available  in  extremely  limited  quantities,  making  direct  measurement  of  its  HSP 
using  a  full  set  of  solvents  impractical.  Although  fluoroalkyl  groups  pose  a  challenge  for  many  group 
contribution  approaches,  general  estimates  of  both  volume  and  cohesive  energy  are  still  possible.  The 
trifluoropropyl  arms  are  also  of  a  small,  prolate  shape,  so  cage  effects  should  be  distinguishable  and  the 
group  contribution  approach  should  be  applicable.  Table  8  lists  the  calculated  volumetric,  cohesive 
energy,  and  HSP  data  for  octakis(trifluoropropyl)  POSS  and  its  components  based  on  the  three 
theoretical  approaches  outlined  earlier,  with  details  of  the  best  available  analog  group  method  provided 
in  Appendix  A.  For  the  Stefanis  approach,  calculation  of  the  cohesive  energy  components  requires  a 
separate  model  for  calculation  of  the  molar  volume.  In  this  case,  the  best  available  analog  group 
method  was  chosen,  since,  in  cases  where  data  on  analogous  compounds  is  readily  available,  it  tended  to 
predict  volumes  more  accurately  (and,  unlike  the  case  for  hexafluoro/.sobutyl  groups,  such  data  is 
readily  available  for  trifluoropropyl  groups).  For  the  HSP  calculations,  the  best  available  analog  group 
method  is  based  on  the  best  available  analog  group  volume  estimates,  the  Beerbower  method  is  based 
on  the  Beerbower  volume  estimates,  and  the  Stefanis  method  does  not  require  a  volume  estimate. 
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Table  8.  Predicted  molar  volume,  cohesive  energy,  and  HSP  for  octakis(trifluoropropyl)  POSS. 


Model  /  Parameter 

Best  Available  Analog  Group 

Beerbower 

Constantinou  /  Stefanis 

Molar  Volume  (cc/mol) 

Each  trifluoropropyl  arm 

79 

79 

n/a 

Cage 

134 

134 

134 

Octakis(trifluoropropyl)  POSS 

766 

766 

n/a 

Cohesive  Energy  (J/mol) 

Dispersive  Component 

Each  trifluoropropyl  arm 

13200 

11300 

19100 

Cage 

65000 

65000 

65000 

Octakis(trifluoropropyl)  POSS 

171000 

156000 

217000 

Polar  Component 

Each  trifluoropropyl  arm 

3940 

4190 

1490 

Cage 

48000 

48000 

48000 

Octakis(trifluoropropyl)  POSS 

80000 

81000 

60000 

Hydrogen  Bonding  Component 

Each  trifluoropropyl  arm 

0 

0 

80 

Cage 

30000 

30000 

30000 

Octakis(trifluoropropyl)  POSS 

30000 

30000 

31000 

Hansen  Solubility  Parameters  (J/cc)12 

Dispersive  Component 

Each  trifluoropropyl  arm 

12.9 

12.0 

17.0 

Cage 

22 

22 

22 

Octakis(trifluoropropyl)  POSS 

14.9 

14.3 

18.2 

Polar  Component 

Each  trifluoropropyl  arm 

7.0 

7.3 

4.8 

Cage 

19 

19 

19 

Octakis(trifluoropropyl)  POSS 

10.2 

10.3 

9.5 

Hydrogen  Bonding  Component 

Each  trifluoropropyl  arm 

0 

0 

1.1 

Cage 

15 

15 

15 

Octakis(trifluoropropyl)  POSS 

6.3 

6.3 

6.9 

As  for  the  radius  of  interaction,  there  is  no  explicit  predictive  method  available.  However,  based  on 
the  data  in  Tables  1  and  2,  it  can  be  seen  that  the  radius  of  interaction  generally  increases  with  the  size 
and  flexibility  of  the  arms.  Also,  comparing  Aobutyl  and  hexafluoro/'.sobutyl  arms,  it  appears  that 
fluorination  leads  to  an  increase  in  the  radius  of  interaction.  The  trifluoropropyl  group  is  somewhat 
smaller  than  the  Aobutyl  group,  but  contains  limited  fluorinated  groups  as  well.  If  these  two  factors 
approximately  counteract  each  other,  one  would  expect  a  radius  of  interaction  for  the 
octakis(trifluoropropyl)  POSS  that  is  similar  to  the  octa(/.sobutyl)  POSS,  or  about  4.3. 

In  order  to  compare  these  predictions  to  experimental  data,  octakis(trifluoropropyl)  POSS  was 
dissolved  in  20  test  solvents.  Using  the  predicted  solubility  parameters  and  the  predicted  radius  of 
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interaction,  RED  values  for  each  solvent  were  assigned  and  compared  against  the  same  “pass/fail” 
criteria  (100  mg/mL  solubility  at  room  temperature)  used  for  the  other  compounds.  The  results  are 
summarized  in  Table  9,  which  shows  both  the  experimental  results  and  codifies  the  outcomes  in  terms 
of  whether  or  not  expectations  were  met  as  well  as  the  nature  of  the  deviations  from  expected  behavior. 


Solvent 

Best  Available  Analog 
Group 

Beerbower 

Stefanis 

RED 

Test 

Code* 

RED 

Test 

Code* 

RED 

Test 

Code* 

1 , 1 ,2,2-tetrachloroethane 

2.18 

Fail 

E- 

2.43 

Fail 

E- 

1.12 

Fail 

E- 

1,4-dioxane 

2.74 

Fail+ 

E- 

2.96 

Fail+ 

E- 

1.83 

Fail+ 

E- 

Acetonitrile 

1.82 

Fail 

E- 

1.85 

Fail 

E- 

2.40 

Fail 

E- 

Acetone 

0.33 

Pass 

E+ 

0.58 

Pass 

E+ 

1.27 

Pass 

U+ 

Allyl  Chloride 

1.64 

Fail 

E- 

1.83 

Fail 

E- 

1.42 

Fail 

E- 

Benzonitrile 

1.38 

Fail+ 

E- 

1.63 

Fail+ 

E- 

0.92 

Fail+ 

U- 

Chloroform 

2.14 

Fail 

E- 

2.33 

Fail 

E- 

1.53 

Fail 

E- 

Cyclohexanone 

1.65 

Fail+ 

E- 

1.90 

Fail+ 

E- 

0.88 

Fail+ 

U- 

Cyclopentane 

2.68 

Fail 

E- 

2.79 

Fail 

E- 

2.64 

Fail 

E- 

Dimethyl  Formamide 

1.83 

Fail 

E- 

2.01 

Fail 

E- 

1.46 

Fail 

E- 

Diphenyl  Ether 

2.66 

Fail 

E- 

2.90 

Fail 

E- 

1.56 

Fail 

E- 

Ethyl  Acetate 

1.23 

Fail+ 

E- 

1.37 

Fail+ 

E- 

1.48 

Fail+ 

E- 

n-Heptane 

2.79 

Fail 

E- 

2.85 

Fail 

E- 

3.05 

Fail 

E- 

Methyl  Ethyl  Ketone 

0.65 

Pass 

E+ 

0.89 

Pass 

E+ 

1.11 

Pass 

U+ 

Nitrobenzene 

2.46 

Fail 

E- 

2.73 

Fail 

E- 

1.08 

Fail 

E- 

N,N-Dimethylacetamide 

1.30 

Pass 

U+ 

1.50 

Pass 

U+ 

1.11 

Pass 

U+ 

o-Dichlorobenzene 

2.30 

Fail 

E- 

2.56 

Fail 

E- 

1.21 

Fail 

E- 

Propylamine 

1.63 

Pass 

U+ 

1.82 

Pass 

U+ 

1.29 

Pass 

U+ 

T  etrahydrofuran 

1.43 

Pass 

U+ 

1.63 

Pass 

U+ 

1.13 

Pass 

U+ 

Toluene 

2.73 

Fail 

E- 

2.92 

Fail 

E- 

2.12 

Fail 

E- 

*a  code  of  “E”  indicates  expected  behavior,  while  “U”  indicates  unexpected  behavior,  the  “+”  indicates  a 
passing  event  and  a  indicates  a  failing  event,  thus  “U-“  for  instance,  indicates  an  unexpected  failure  (i.e.  lack 
of  expected  solubility). 


Examining  Table  9,  several  factors  should  be  kept  in  mind.  First,  the  group  contributions  for  the 
POSS  cage  were  derived  from  the  data  on  aromatic  rather  than  aliphatic  or  fluorinated  POSS 
compounds,  yet  they  appear  to  have  predicted  the  properties  of  the  fluoroaliphatic  POSS  compound 
quite  well.  Second,  the  differences  in  the  performance  of  the  three  methods  (with  the  best  available 
analog  group  method  and  the  Beerbower  method  appearing  equally  proficient)  stem  primarily  from 
different  predictions  of  the  dispersive  component  of  the  fluoroalkyl  arm  (see  Table  8),  with  the  Stefanis 
method  being  the  highest  estimate.  Therefore,  from  a  group  contribution  standpoint,  the  presence  of 
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POSS  cages  appears  not  to  be  a  limiting  factor.  Some  improvement  in  predictive  capability  in  the  case 
of  the  Stefanis  method  can  be  obtained  by  using  a  larger  radius  of  interaction. 

Based  on  the  limited  data  set  in  Table  9,  a  maximization  of  %  according  to  the  method  described  in 

1/9  1/9 

the  Experimental  section  yields  centroid  parameters  of  8d  =  16.9  (J/cc)  ,  5p  =  9.1  (J/cc)  ,  8h  =  8.9 
(J/cc)  ,  and  Ro  =  4.5  with  radii  of  gyration  of  0.3  -  0.5  (J/cc)  ,  five  good  solvents  and  no  exceptions  to 
the  “spherical  rule”.  Given  the  limited  data  set  and  the  many  uncertainties  inherent  in  group 
contribution  methods,  the  predicted  values  (which  involved  no  adjustable  parameters)  agreed  quite  well 
with  the  experimental  data.  These  results  reinforce  the  earlier  results  that  suggest  the  HSP  approach  is 
well-suited  to  POSS  cages,  and  that  group  contributions  for  POSS  cages  can  be  utilized  on  at  least 
certain  types  of  POSS.  They  also  illustrate  that  the  quality  of  available  group  contribution  values  for  the 
peripheral  arms  on  POSS  cages  represents  a  major  consideration  when  developing  group  contributions 
for  the  cage. 

Although  the  present  work  has  shown  the  utility  of  the  HSP  approach  in  understanding  the  solubility 
behavior  of  POSS  compounds  in  small  molecule  solvents,  a  significantly  greater  payoff  will  be  achieved 
if  these  results  can  be  extended  to  POSS/polymer  blends.  In  considering  the  HSP  approach  to 
predicting  the  miscibility  of  POSS  compounds  with  polymers,  there  remain  several  important,  yet 
unanswered  questions.  First,  in  the  present  work,  as  when  determining  HSP  for  polymers,  the  radius  of 
interaction  is  considered  a  characteristic  of  the  POSS  compound  rather  than  the  solvent.  For 
POSS/polymer  blends,  however,  it  is  unclear  which  radius  is  more  appropriate  -  that  of  the  polymer, 
that  of  the  POSS  compound  (which  does  vary  significantly),  or  some  combination  of  both.  Second,  a 
widely  reported  benefit  to  blending  POSS  compounds  into  polymers  is  the  creation  of  free  volume, 
which  presumably  results  from  an  inability  of  polymer  chains  to  penetrate  the  complex  geometry  of  the 
POSS  nanostructure.  In  polymer/solvent  and  POSS/solvent  systems,  the  relatively  small  size  and 
comparatively  simple  geometry  of  the  solvent  molecules  allows  them  to  “fit  in”  to  complex  geometries 
and  mitigate  these  effects,  but  in  POSS/polymer  systems,  the  “intermeshing”  of  nanostructures,  or  lack 
thereof,  may  play  a  significant  role  in  controlling  miscibility.  The  HSP  approach  is  not  well-suited  to 


26 


capturing  these  potentially  important  effects.  The  present  work,  however,  by  demonstrating  that  HSP 
values  are  a  meaningful  representation  of  cohesive  energy  in  POSS  compounds  (and  potentially  for 
silicate  nanostructures  in  general),  are  compatible  with  group  contribution  approaches,  and  have 
significant  predictive  power,  effectively  lays  the  groundwork  for  future  investigations  of  these 
remaining  issues. 

4.  Conclusions 

Hansen  Solubility  Parameters  (HSP)  for  five  different  Polyhedral  Oligomeric  SilSesquioxane  (POSS) 
compounds  have  been  successfully  determined,  demonstrating  that  the  HSP  approach  developed 
primarily  for  organic  compounds  and  polymers  may  be  extended  to  at  least  limited  types  of  organic- 
inorganic  compounds.  In  general,  the  relative  energy  difference  (i.e.  “spherical  rule”)  provided  good  to 
excellent  predictive  power,  with  somewhat  less  reliable  predictions  for  octa(Aobutyl)  POSS  at  100 
mg/mL.  Though  a  consistent  analysis  procedure  was  employed  for  all  compounds,  it  was  found  that 
alternative  analyses  provided  superior  estimates  of  HSP  when  either  the  number  of  good  solvents  was 
very  low,  or  when  too  few  poor  solvents  were  available  to  properly  constrain  the  boundary  of  the  sphere 
of  solubility.  Group  contributions  for  the  octameric  POSS  cage  were  determined  using  three  different 
approaches,  which  gave  similar  results,  with  the  best  estimate  being  5d  =  22  (J/cc)  ,  8P  =  19  (J/cc)  , 
and  5h  =  15  (J/cc)  ,  with  an  estimated  uncertainty  of  approx.  5  (J/cc)  .  A  major  constraint  in  deriving 
and  using  these  estimates  was  the  ability  of  the  various  group  contribution  methods  to  correctly  model 
the  contributions  from  the  peripheral  arm  groups.  The  utility  of  the  HSP  approach  was  demonstrated  by 
successfully  identifying  mixtures  of  poor  solvents  that  provided  significantly  enhanced  solubility  for 
octa(Aobutyl)  POSS,  and  by  successfully  estimating  the  HSP  of  octakis(trifluoropropyl)  POSS  from 
group  contributions  derived  solely  from  aromatic  POSS  compounds. 
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Appendix  A:  Description  of  the  Best  Available  Analog  Group  Method 

Al.  Estimates  of  Molar  Volume 

Although  many  theoretical  methods  allow  for  direct  estimation  of  the  cohesive  energy  of  a  given 
chemical  substance,  certain  methods,  such  as  the  group  contribution  approach  of  Stefanis,  involve 
direct  calculation  of  solubility  parameters.  These  methods  often  do  not  produce  results  that  are 
independent  of  molecular  size.  For  small  molecules,  these  size  effects  must  be  considered,  but  for  very 
large  molecules,  such  approaches  provide  divergent  results.  A  method  of  appropriately  combining  the 
results  for  smaller  fragments  must  be  developed.  Because  the  cohesive  energy  is  an  extensive  variable, 
it  forms  the  natural  basis  for  additive  combination.  Therefore,  for  large  molecules,  a  more  suitable 
approach  is  to  convert  the  calculated  solubility  parameters  for  appropriately  sized  fragments  into 
cohesive  energies,  add  these  together,  and  then  calculate  the  solubility  parameter  for  the  large  molecule 
as  a  whole.  Such  an  approach  was  followed  for  POSS  compounds. 

The  aforementioned  approach  requires  a  means  of  estimating  the  molar  volume  of  molecular 
fragments.  Although  various  correlations,  such  as  those  from  Beerbower  as  reported  by  Hansen1  and 
from  Constantinou,65  may  be  used  for  this  purpose,  they  are  often  limited  to  only  the  most  common 
organic  chemical  functionalities.  In  particular,  as  discussed  in  the  main  body  of  the  paper,  the  molar 
volumes  of  naphthyl  and  fluoroalkyl  groups  found  on  the  POSS  compounds  under  study  are  not 
calculable  using  one  or  more  of  the  above  methods.  In  order  to  provide  another  basis  for  calculation, 
the  “best  available  analog  group”  method  was  adopted.  Although  this  method  lacks  a  consistent  set  of 
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defined  rules,  it  is  based  on  straightforward  principles,  and  can  be  applied  to  almost  any  chemical 
structure.  It  also  has  the  significant  advantage  of  being  based  on  existing  data  and/or  calculation. 

Put  simply,  the  “best  available  analog  group”  method  involves  finding  the  closest  chemical  analog  of 
the  group  of  interest  for  which  either  experimental  data  or  more  well-established  estimates  are  already 
available,  and  using  general  rules  to  translate  these  values  into  new  estimates.  For  instance,  for  the 
phenethyl  group  found  in  octa(phenethyl)  POSS,  the  compound  ethylbenzene  is  a  close  chemical  analog 
with  experimentally  determined  volume  (and  cohesive  energy)  values.  One  needs  only  to  apply  a 
general  rule  for  translating  the  molar  volume  of  ethylbenzene  into  the  molar  volume  of  the  phenethyl 
group  (an  ethylbenzene  with  the  terminal  aliphatic  hydrogen  replaced  by  a  chemical  bond  to  a  silicon  of 
the  POSS  cage). 

The  effect  of  substituting  a  bond  to  a  POSS  silicon  for  a  terminal  hydrogen  in  general  is  not  readily 
available.  However,  a  simpler  rule  can  be  used  to  provide  a  rough  estimate  based  on  the  properties  of 
dimeric  compounds.  For  instance,  the  dimer  of  benzene  is  biphenyl.  Biphenyl  may  be  visualized  as 
consisting  of  two  phenyl  groups  joined  to  one  another,  thus  one  means  of  estimating  the  molar  volume 
of  a  phenyl  group  is  simply  to  take  half  the  molar  volume  of  biphenyl.  For  many  aliphatic  compounds, 
the  molar  volume  of  the  dimer  is  about  35  cc/mol  less  than  twice  the  volume  of  the  monomer,  that  is  to 
say,  the  shrinkage  on  dimerization  is  about  35  cc/mol.  (Numerous  examples  are  provided  in  Table  S4  in 
Supporting  Information.)  For  a  group  like  phenethyl,  with  the  bonded  terminal  on  an  aliphatic  group, 
the  molar  volume  may  be  estimated  as  17.5  cc/mol  less  than  the  volume  of  the  corresponding  non- 
bonded  compound.  For  the  phenethyl  example,  the  molar  volume  estimate  would  be  17.5  cc/mol  less 
than  the  molar  volume  of  ethylbenzene.  Note  that  for  aromatic  linkages,  the  shrinkage  on  dimerization 
is  somewhat  lower,  at  about  12  cc/mol. 

To  calculate  the  molar  volume  of  a  POSS  compound  by  the  “best  available  analog  group”  method, 
one  simply  needs  to  identify  the  analog  group,  determine  the  shrinkage  on  dimerization,  and  then 
subtract  half  this  value  from  the  molar  volume  of  the  analog  group  to  estimate  the  volume  of  each  arm. 
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The  volume  of  the  cage  itself  may  then  be  added  based  on  the  model  geometry  estimate  described  in 
Appendix  B.  For  all  the  POSS  cage  types  studied,  the  shrinkage  on  dimerization  was  based  on  the 
presence  of  either  an  aliphatic  or  aromatic  linkage  at  the  root  of  each  arm.  For  phenethyl,  styrenyl, 
phenyl,  1-napthyl,  and  /.sobutyl  arms,  the  molar  volumes  of  the  analogs  (ethylbenzene,  styrene,  benzene, 
naphthalene,  and  /.sobutanc)  were  all  readily  available.  For  the  trifluoropropyl  ann,  the  volume  of  the 
analog  (1,1,1-trifluoropropane)  is  also  available.  For  both  /.sobutanc  and  1,1,1-trifluoropropane,  the 
volume  of  the  saturated  liquid  at  298  K  (which  is  at  higher  than  atmospheric  pressure)  was  used.  Note 
that  all  estimates  and  data  sources  are  tabulated  in  Supporting  Information. 

For  the  hcxafluoro/.vobutyl  arm,  reliable  volumetric  data  on  the  analog  compound  1, 1,1, 3,3,3- 
hexafluoro-2-methylpropane  was  not  readily  available.  However,  volumetric  data  for  the  unsaturated 
analog  (hexafluoro/sobutylene)  was  readily  available.  In  this  case,  another  general  rule  needed  to  be 
applied,  based  on  the  expected  effect  of  hydrogenation.  The  expected  effect  of  hydrogenation  was 
determined  by  subtracting  the  molar  volume  of  isobutylene  from  the  molar  volume  of  /'.sobutanc, 
resulting  in  an  estimate  of  a  9  cc/mol  increase  in  volume  due  to  hydrogenation.  This  value  was  added  to 
the  experimental  molar  volume  of  hcxafluoro/'.sobutylcnc  to  obtain  the  estimated  molar  volume  of 
1,1,1 ,3,3,3-hexafluoro/.sobutanc,  which  was  then  converted  to  the  estimated  volume  of  the 
corresponding  ann  by  subtracting  half  the  shrinkage  on  dimerization  for  an  aliphatic  group  (17.5 
cc/mol)  to  obtain  the  estimated  arm  volume.  Thus,  the  estimated  volume  for  the  hexafluoro/.vobutyl  arm 
is  less  reliable  than  for  the  arms  with  more  common  analogs,  and  involves  certain  arbitrary  choices, 
illustrating  both  the  strength  (basis  in  experimental  data)  and  weakness  (lack  of  a  straightforward 
universal  step-by-step  formula  for  computation)  inherent  in  this  approach. 

A2.  Estimates  of  Cohesive  Energy  Components 

Given  the  complexities  involved  in  calculating  a  quantity  as  simple  as  molar  volumes  by  the  “best 
available  analog  group”  method,  it  might  seem  at  first  that  calculation  of  cohesive  energy  components 
would  be  even  more  challenging.  On  the  contrary,  such  calculations  are  in  fact  much  simpler,  because 
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the  foundation  of  many  successful  group  contribution  methods  for  estimation  of  the  molar  cohesive 
energy  is  the  successful  assumption  that  it  is  directly  proportional  to  the  number  of  functional  groups 
with  only  a  small  dependence  on  how  those  groups  are  arranged.  For  instance,  the  cohesive  energy  of 
the  compound  biphenyl  is  assumed  to  be  exactly  twice  that  of  benzene  in  numerous  group  contribution 
methods  that  have  proven  reliable  over  many  years.  Therefore,  attachment  to  a  POSS  cage  can  be 
assumed  to  have  only  a  small  effect  on  the  cohesive  energy  of  a  particular  arm,  and  thus  the  cohesive 
energy  of  a  given  arm  may  be  assumed  to  equal  that  of  its  free  molecule  analog  (e.g.  the  cohesive 
energy  of  a  phenethyl  arm  is  assumed  to  equal  that  of  ethylbenzene).  Moreover,  for  the  purposes  of  the 
method,  it  is  assumed  that  the  distribution  of  cohesive  energy  components  also  does  not  change 
significantly  with  attachment  to  a  POSS  cage. 

For  phenethyl,  styrenyl,  naphthyl,  and  phenyl  arms,  experimentally  derived  cohesive  energy 
components  are  widely  available,  and  for  computational  purposes,  the  values  provided  by  Hansen  were 
used.  For  /.vobutyl  anns,  the  cohesive  energy  components  of  /.vobutanc  were  needed.  In  general, 
aliphatic  hydrocarbons  are  assumed  to  have  no  polar  or  hydrogen  bonding  components  of  cohesive 
energy.  Therefore,  for  Aobutane,  the  total  cohesive  energy  was  determined  by  the  thermodynamic 
relation 


E  =  AHvap  -  RT  (A-l) 

where  AHvap  is  the  molar  heat  of  vaporization  at  298  K,  the  experimental  value  of  which  was  readily 
available,  R  is  the  universal  gas  constant,  and  T  the  absolute  temperature.  The  total  cohesive  energy 
was  then  assumed  to  equal  the  dispersive  component  of  the  cohesive  energy,  with  the  other  components 
being  zero.  Note  that  all  estimates  and  data  sources  are  tabulated  in  Supporting  Information. 

For  the  hexafluoroAobutyl  arm,  as  mentioned  previously,  experimental  data  on  the  analogous 
compound,  hexafluoroAobutylene,  was  not  readily  available.  However,  the  components  of  cohesive 
energy  for  hexafluoroAobutylene,  were  estimated  in  the  tables  provided  by  Hansen.1  As  when 
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estimating  molar  volume,  comparative  heat  of  vaporization  data  for  isobutylene  and  /.sobutanc  were  also 
used  to  estimate  the  effect  of  hydrogenation  on  the  dispersive  component  of  the  cohesive  energy,  and 
the  difference  was  applied  as  a  correction  to  the  value  for  hexafluoro/'.sobutylcne  to  estimate  the  value 
for  h cx a (1  uo ro /.sob utan e .  Because  the  dispersive  energy  of  organic  compounds  is  often  calculated  on  the 
basis  of  the  total  cohesive  energy  of  the  hydrocarbon  homo-morph  for  non-fluorinated  compounds,  and 
because  the  difference  in  saturation  involves  a  non-fluorinated  portion  of  the  molecule,  the 
aforementioned  procedure  appears  reasonable  for  the  dispersive  component  of  the  cohesive  energy.  For 
the  polar  and  hydrogen  bonding  components,  a  quick  scan  of  cohesive  energy  data  tables  indicated  that 
no  obvious  simple  correction  exists  for  comparing  saturated  and  unsaturated  compounds,  so  no 
estimation  of  these  components  was  attempted  for  the  hexafluoro/sobutyl  arm.  Rather,  they  were 
simply  assumed  to  equal  zero.  Note  that  the  hydrogen  bonding  component  of  the  cohesive  energy  for 
many  fluorocarbons  is  generally  assumed  to  equal  zero  in  many  group  contribution  methods  as  well, 
leaving  only  the  polar  component  in  question.  Because  fluorination  of  aliphatic  hydrocarbons  affects 
both  dispersive  and  polar  components  of  cohesive  energy,  the  hydrocarbon  homo-morph  does  not 
provide  a  reliable  means  of  estimating  only  the  dispersive  component. 

For  the  trifluoropropyl  arm,  cohesive  energy  components  for  1,1,1-trifluoropropane  were  not  readily 
available,  though  an  experimental  value  for  AHvap  at  298  K  was  found.  The  closest  analog  with 
component  data  estimated  by  Hansen  is  1,1,1-trifluoroethane,  although  based  on  the  NIST 
thermodynamic  tables  for  1,1,1-trifluoroethane,  it  appears  the  molar  volume  and  solubility  parameters 
listed  are  valid  only  at  174  K.  The  same  NIST  thermodynamic  tables,  however,  allow  for  the 
calculation  of  AHvap  at  298  K.  Although  the  distribution  of  cohesive  energy  among  its  components 
could  vary  somewhat  with  temperature,  for  simplicity  it  was  assumed  that  the  distribution  remained 
unchanged,  and  the  component  values  were  simply  scaled  to  the  corrected  total.  Having  estimated 
cohesive  energy  components  for  1,1,1-trifluoroethane,  it  was  assumed  that  the  presence  of  an  additional 
methylene  group  would  add  mainly  to  the  dispersive  component  of  the  cohesive  energy.  Therefore,  the 
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polar  and  hydrogen  bonding  components  were  taken  to  be  those  of  1,1,1-trifluoroethane,  while  the 
difference  between  their  sum  and  the  experimentally  determined  total  based  on  AHvap  at  298  K  was 
assigned  to  the  dispersive  component. 

These  examples  illustrate,  once  again,  the  need  for  judicious  intuition  that  makes  the  best  available 
analog  group  method  less  attractive  in  some  respects.  In  addition,  it  should  be  noted  that  in  some  cases 
the  values  for  the  analog  compounds  themselves  are  merely  better  established  estimates.  Therefore, 
while  the  best  available  analog  group  method  has  the  advantage  of  a  more  direct  connection  to 
experimental  data  for  common  compounds,  in  cases  such  as  fluorocarbons,  where  other  group 
contribution  approaches  have  limited  success,  this  method  also  suffers  from  significant  limitations. 

Appendix  B:  Estimate  of  the  Molar  Volume  of  a  POSS  Cage  Based  on  Molecular  Geometry 

Although  it  is  possible  to  estimate  the  molar  volume  contribution  of  a  POSS  cage  by  comparing 
estimates  for  the  arm  volume  with  experimentally-derived  total  volume  values,  it  is  useful  to  have  a 
theoretical  comparison  because  the  theoretical  arm  values  contain  numerous  uncertain  assumptions, 
while  the  experimentally-derived  values  are  for  powders  of  often  unknown  extent  of  crystallinity.  To 
determine  a  theoretical  value,  the  molecular  geometry  of  the  octameric  POSS  cage  may  be  utilized. 
Although  the  exact  details  of  the  cage  geometry  differ  among  various  compounds  and  are  likely  to 
change  slightly  based  on  the  ordering  of  their  surroundings,  for  the  purposes  outlined  above,  they  are  of 
sufficient  similarity  to  be  considered  as  a  single  entity. 

In  order  to  construct  a  geometric  model,  three  forms  of  information  are  needed.  The  first  is  the 
relative  arrangement  of  the  centers-of-mass  of  each  atom  in  the  cage.  Such  data  is  readily  available  for 
X-ray  studies  of  POSS  crystals,66  and  the  data  for  octa(methyl)  POSS  were  chosen  as  a  representative 
example.  The  second  is  the  size  of  atoms.  For  this  purpose,  the  van  der  Waals  radii,  as  commonly 
reported,  were  chosen.  The  third,  and  most  challenging,  is  an  appropriate  bounding  shape  that  sets  apart 
the  region  of  space  that  “belongs  to”  the  POSS  cage. 
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There  are  many  possible  approaches  to  determining  an  appropriate  boundary.  For  instance,  many- 
atom  molecular  models  of  POSS  compounds  in  various  liquids  could  be  constructed.  From  these,  a 
Voronoi  tessellation  using  each  atom  could  be  constructed,  and  the  volume  of  the  Voronoi  polygons 
assigned  to  atoms  of  the  POSS  cage  totaled  and  averaged  over  many  models.  Such  an  approach, 
though,  is  highly  complex  for  the  aforementioned  purpose  of  simply  providing  a  reasonable  estimate  to 
compare  with  experimental  data.  A  much  simpler  approach  is  simply  to  consider  the  pertinent 
geometric  factors  and  principles  that  would  determine  an  effective  boundary.  The  latter,  simpler 
approach  was  chosen. 

Although  POSS  compounds  are  known  to  create  additional  free  volume  in  polymers  and  liquids,  in 
general  the  packing  fraction  in  these  liquids  is  still  expected  to  be  close  to  the  maximum.  Therefore, 
any  boundary  should  closely  conform  to  the  van  der  Waals  radii  of  the  cage.  Although  the  cage  itself 
does  have  certain  slightly  concave  features,  in  a  random  arrangement  of  molecules,  not  many  are  likely 
to  fit  snugly  into  such  features,  and  a  reasonable  assumption  is  that  any  boundary  should  define  a 
convex  interior  region.  Because  the  octameric  POSS  cage  is  approximately  cubic,  the  simplest 
boundary  would  be  a  cube  that  circumscribed  the  atoms  in  the  cage  (as  defined  by  their  van  der  Waals 
radii),  which  equates  to  having  an  edge  length  of  0.684  mn.  Figure  B1  shows  such  a  cube  superimposed 
in  a  space-filling  model  of  the  octameric  POSS  cage. 

The  main  difficulty  with  such  a  cube  is  that  the  center  of  mass  of  the  peripheral  carbon  atom  nearest  to 
each  silicon  atom  actually  falls  within  the  cube,  so  each  corner  of  the  cube  contains  a  region  that  is  more 
properly  associated  with  the  peripheral  groups.  To  correct  this  issue,  the  cube  can  be  truncated  along  a 
plane  normal  to  the  [111]  direction  (using  the  principal  axes  of  the  cube  as  principal  directions). 
Although  the  plane  that  bisects  the  silicon-carbon  bond  seems  a  natural  choice,  one  must  consider  that 
the  bond  length  in  question  varies,  and  that  the  carbon  atom  is  smaller  than  the  silicon  atom,  such  that 
the  overlap  of  spheres  defined  by  van  der  Waals  radii  is  not  centered  on  the  bisecting  plane.  An 
“internally  consistent”  choice  for  the  octameric  POSS  cage  would  be  the  plane  that  is  defined  by  the 
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three  points  near  each  vertex  of  the  cube  where  the  van  der  Waals  radius  of  the  silicon  atom  makes 
contact  with  the  cube.  Although  this  plane  is  likely  to  pass  near  the  center  of  mass  of  each  of  the  nearest 
peripheral  carbon  atoms,  and  thus  may  attribute  too  much  volume  to  the  cage  near  the  comers,  a  more 
accurate  solution  would  require  a  concave  shape  in  order  not  to  exclude  an  inordinate  amount  of  volume 
within  the  van  der  Waals  radius  of  the  silicon  atoms.  The  location  of  the  plane  used  to  truncate  the  cube 
is  also  shown  in  Figure  B 1 . 

When  the  appropriate  truncation  is  applied  to  each  comer  of  the  cube,  the  resultant  multi-faceted 
shape  has  an  appearance  given  by  the  rough  sketches  in  Figure  B2.  The  region  is  characterized  by 
diamond-shaped  faces  coincident  with  the  cube  that  circumscribes  the  cage,  and  by  six-sided  truncated 
triangular  faces  formed  by  the  planes  that  tmncate  each  corner  of  the  cube.  Figure  B2  also  includes  a 
superimposed  chemical  structure  of  a  portion  of  the  cage  to  help  clarify  the  relationships  between  the 
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geometry  of  the  faces  and  positions  of  the  atoms  in  the  cage.  The  total  volume  of  the  cube  (0.3200  nm  ) 
equates  to  192.7  cc/mol,  while  that  of  the  regions  excluded  by  truncation  (0.0969  nm  ),  or  58.4  cc/mol, 
which  leaves  134  cc/mol  as  the  volume  assigned  to  the  cage. 

To  check  this  value,  the  difference  between  the  experimental  volume  of  the  POSS  compounds  studied 
and  the  estimated  arm  volume  (obtained  by  averaging  the  results  of  the  three  models  shown  in  Table  3) 
can  be  calculated.  The  results  are:  for  octa(methyl)  POSS,  129  cc/mol,  for  octa(phenethyl)  POSS  139 
cc/mol,  for  octa(styrenyl)  POSS,  166  cc/mol,  for  octa(Aobutyl)  POSS  123  cc/mol,  for  (1-naphthyl)- 
heptaphenyl  POSS,  188  cc/mol,  and  for  octakis(hexafhroroAobutyl)  POSS  (for  which  the  models  are 
less  reliable),  -20  cc/mol.  Excluding  the  latter,  non-physical  result,  the  average  is  149  cc/mol,  with  a 
standard  deviation  of  27  cc/mol.  Thus,  the  model  does  reasonably  well  in  describing  the  experimental 
data.  Furthermore,  of  the  various  experimental  volumes  calculated,  the  octa(methyl)  POSS  is  likely  to 
involve  the  least  amount  of  relative  error  because  it  represents  the  smallest  peripheral  group 
contribution,  and  provides  the  most  reliable  estimate  of  the  cage  volume.  From  that  standpoint,  the 
model  also  provides  a  reasonable  estimate. 
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0.684  nm 


Figure  Bl.  Key  features  of  the  geometric  model  of  an  octameric  POSS  cage,  including 
circumscribing  cube  and  the  plane  used  to  truncate  the  corners  in  order  to  account  for  the  presence  of 
peripheral  groups. 


Figure  B2.  Sketches  of  the  truncate  cube  used  to  define  the  effective  volume  contribution  from  a 
octameric  POSS  cage,  including  (left)  view  along  the  principal  axes  of  the  cube  shown  in  Figure  Bl, 
and  (right)  rough  sketch  of  the  isometric  view  showing  the  chemical  structure  of  a  portion  of  the  cage 
for  reference. 
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Table  SI  lists  the  names  of  all  45  solvents  used  in  the  “standard”  test  set  for 
determination  of  the  Hansen  Solubility  Parameters  (HSP)  of  five  octameric  POSS  compounds, 
octa(phenethyl),  octa(styrenyl),  octa(Aobutyl),  octa(hexafluoroAobutyl),  and  (1-naphthyl)- 
heptaphenyl  POSS.  Note  that  due  to  limited  supplied  of  material,  only  36  of  the  45  solvents 
were  used  with  octa(hexafluoroAobutyl)  POSS.  Table  1  provides  both  the  test  results  and 
whether  or  not  those  results  agree  with  the  “spherical  rule”  based  on  the  subsequent  HSP 
determination.  Details  of  the  experiments,  the  ratings,  and  the  HSP  determination  procedure 
may  be  found  in  the  main  body  of  the  paper.  “Pass”  indicates  total  dissolution  within  two 
minutes,  “Pass  — “  indicates  total  dissolution  within  two  minutes  to  one  hour,  “Fail”  indicates  no 
dissolution  whatsoever,  and  “Fail+”  indicates  incomplete  dissolution  with  some  visible  signs  of 
good  dispersion  or  partial  dissolution.  Note  that  results  of  the  “Supplemental”  testing  are  listed 
in  Table  S2  (see  the  main  body  of  the  text  for  a  detailed  description  of  supplemental  tests). 
Briefly,  whereas  the  “standard”  tests  were  carried  out  at  100  mg/mL,  test  PE-A1  utilized  200 
mg/mL,  and  test  NP-A1  utilized  50  mg/mL.  Test  PE-A2  utilized  additional  solvents  and  solvent 
mixtures,  listed  in  Table  S3. 


Approved  for  public  release,  distribution  unlimited. 


Table  SI 


Solubility  Test  Results  for  “Standard”  Tests 


Solvent 

PE 

ST 

IB 

HF 

NP 

P/F 

E/U 

P/F 

E/U 

P/F 

E/U 

P/F 

E/U 

P/F 

E/U 

1,1,2,2-Tetrabromoethane 

P- 

E 

F+ 

U 

F+ 

E 

F 

E 

F 

E 

1,1,2,2-Tetrachloroethane 

P 

E 

P 

E 

P- 

E 

F 

E 

F+ 

E 

1,4-Butanediol 

F 

E 

F 

E 

F 

E 

F 

E 

F 

E 

1,4-Dioxane 

P- 

E 

F 

U 

F+ 

E 

F 

E 

F 

E 

1-Decene 

F+ 

E 

F 

E 

F 

E 

F 

E 

2-Chlorophenol 

P 

E 

P 

E 

F 

E 

F 

E 

2-Ethyl- 1  -hexanol 

F+ 

E 

F 

E 

F 

E 

F 

E 

Acetonitrile 

F 

E 

F 

E 

F 

E 

F 

E 

F 

E 

Acetone 

P 

E 

F+ 

E 

F+ 

E 

P 

E 

F 

E 

Allyl  chloride 

P 

E 

P 

E 

F 

E 

F 

E 

F 

E 

Benzene 

P 

E 

P 

E 

F+ 

U 

F 

E 

Benzonitrile 

P 

E 

P- 

E 

P 

E 

F 

E 

F+ 

E 

Benzyl  chloride 

P 

E 

P 

E 

P 

E 

F 

E 

F 

E 

Chloroform 

P 

E 

P 

E 

P 

E 

F 

E 

P 

E 

Cyclohexane 

P- 

U 

F 

E 

P 

E 

F 

E 

F 

E 

Cyclohexanone 

P 

E 

F+ 

U 

F+ 

E 

F 

E 

F+ 

E 

Cyclopentane 

F+ 

E 

F 

E 

P 

E 

F 

E 

F 

E 

Diethyl  ether 

P 

U 

F 

E 

P 

U 

F 

E 

F 

E 

Dimethyl  Formamide 

P- 

E 

P 

E 

F 

E 

P 

E 

F 

E 

Dimethyl  Sulfoxide 

F 

E 

F 

E 

F 

E 

F 

E 

F 

E 

Di-n-butyl  ether 

F+ 

E 

F 

E 

F 

E 

F 

E 

Diphenyl  Ether 

P- 

E 

P- 

E 

F+ 

E 

F 

E 

F+ 

E 

Dodecane 

F+ 

E 

F 

E 

F+ 

E 

F 

E 

Ethanol 

F 

E 

F 

E 

F 

E 

F 

E 

F 

E 

Ethyl  Acetate 

P 

E 

P 

E 

F 

E 

P 

E 

F 

E 

EGDE 

P 

E 

P 

E 

F 

E 

P 

E 

F 

E 

n-Heptane 

F 

E 

F 

E 

F 

E 

P 

U 

F 

E 

Hexadecane 

F 

E 

F 

E 

F+ 

E 

F 

E 

Isopropanol 

F 

E 

F 

E 

F 

E 

F 

U 

F 

E 

Methyl  Ethyl  Ketone 

P 

E 

P 

E 

F 

E 

P 

E 

F 

E 

n-Butanol 

F 

E 

F 

E 

F 

E 

P 

E 

F 

E 

n-Hexane 

F 

E 

F 

E 

P 

U 

F 

E 

F 

E 

Nitrobenzene 

P 

E 

P 

E 

F+ 

E 

F 

E 

F 

E 

N-methyl-2-pyrrolidone 

P 

E 

P 

E 

F+ 

E 

P 

U 

N-N-dimethylacetamide 

P 

E 

F+ 

U 

F+ 

E 

P 

E 

F 

E 

n-Propanol 

F 

E 

F 

E 

F 

E 

P 

E 

F 

E 

o-Dibromobenzene 

P 

E 

P- 

E 

F 

E 

F 

E 

o-Dichlorobenzene 

P 

E 

F 

U 

F+ 

E 

F 

E 

F 

E 

Phenyl  acetylene 

P 

E 

P 

E 

F 

U 

F 

E 

F 

E 

Propylamine 

P 

E 

P 

E 

F 

E 

P 

E 

F 

E 

2 


Solvent 

PE 

ST 

IB 

HF 

NP 

P/F 

E/U 

P/F 

E/U 

P/F 

E/U 

P/F 

E/U 

P/F 

E/U 

Propylene  Carbonate 

F+ 

E 

F 

E 

F 

E 

F 

E 

F 

E 

Quinoline 

P- 

E 

P 

E 

F 

E 

F 

E 

F 

E 

Tetrahydrofuran 

P 

E 

P 

E 

P 

U 

P 

E 

P 

E 

Toluene 

P 

E 

P 

E 

P- 

E 

F 

E 

F 

E 

Xylenes 

P 

E 

P- 

E 

F+ 

U 

F 

E 

F 

E 

Guide  to  abbreviations: 

EGDE:  ethylene  glycol  dimethyl  ether 

PE:  octa(phenethyl)  POSS 

ST:  octa(styrenyl)  POSS 

IB:  octa(z'sobutyl)  POSS 

HF:  octak i s(h cxa  11  uoro/.vobutyl )  POSS 

NP:  (l-naphthyl)heptaphenyl  POSS 

E:  as  expected,  based  on  best  fit  HSP  and  “spherical  rule” 

U:  not  as  expected,  based  on  best  fit  HSP  and  “spherical  rule” 
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Table  S2 


Solubility  Test  Results  for  “Supplemental”  Tests 


Solvent 

PE-A1 

NP-A1 

NP-A2 

P/F 

E/U 

P/F 

E/U 

P/F 

E/U 

1,1,2,2-Tetrabromoethane  P 

E 

F 

E 

F 

E 

1,1,2,2-Tetrachloroethane  P 

E 

F+ 

E 

F+* 

E 

1,4-Butanediol 

F 

E 

F 

E 

F 

E 

1,4-Dioxane 

P 

E 

F 

E 

F 

E 

1-Decene 

F+ 

E 

F 

E 

F 

E 

2-Chlorophenol 

F 

U 

F 

E 

F 

E 

2-Ethyl- 1  -hexanol 

F+ 

E 

F 

E 

F 

E 

Acetonitrile 

F 

E 

F 

E 

F 

E 

Acetone 

F 

E 

F 

E 

F 

E 

Allyl  chloride 

F 

E 

F 

E 

F 

E 

Benzene 

F+ 

E 

F 

E 

F 

E 

Benzonitrile 

P 

E 

F+ 

E 

F+* 

E 

Benzyl  chloride 

P 

E 

F 

E 

F 

E 

Chloroform 

P 

E 

P 

E 

P 

U 

Cyclohexane 

F 

E 

F 

E 

F 

E 

Cyclohexanone 

P 

E 

F+ 

E 

F+* 

E 

Cyclopentane 

F 

E 

F 

E 

F 

E 

Diethyl  ether 

F+ 

E 

F 

E 

F 

E 

Dimethyl  Formamide 

P 

E 

F 

E 

F 

E 

Dimethyl  Sulfoxide 

F 

E 

F 

E 

F 

E 

Di-n-butyl  ether 

F+ 

E 

F 

E 

F 

E 

Diphenyl  Ether 

F 

U 

F+ 

E 

F+* 

U 

Dodecane 

F+ 

E 

F 

E 

F 

E 

Ethanol 

F 

E 

F 

E 

F 

E 

Ethyl  Acetate 

F 

E 

F 

E 

F 

E 

EGDE 

F 

E 

F 

E 

F 

E 

n-Heptane 

F 

E 

F 

E 

F 

E 

Hexadecane 

F 

E 

F 

E 

F 

E 

Isopropanol 

F 

E 

F 

E 

F 

E 

Methyl  Ethyl  Ketone 

F 

E 

F 

E 

F 

E 

n-Butanol 

F 

E 

F 

E 

F 

E 

n-Hexane 

F 

E 

F 

E 

F 

E 

Nitrobenzene 

P- 

E 

F 

E 

F 

E 

N-methyl-2-pyrrolidone 

P- 

E 

P 

U 

P 

E 

N-N-dimethylacetamide 

P- 

E 

F 

E 

F 

E 

n-Propanol 

F 

E 

F 

E 

F 

E 

o-Dibromobenzene 

P 

E 

F 

E 

F 

E 

o-Dichlorobenzene 

P 

E 

F 

E 

F 

E 

Phenyl  acetylene 

P 

E 

F 

E 

F 

E 

Propylamine 

P 

E 

F 

E 

F 

E 

Solvent 

PE-A1 

NP-A1 

NP-A2 

P/F 

E/U 

P/F 

E/U 

P/F 

E/U 

Propylene  Carbonate 

F 

E 

F 

E 

F 

E 

Quinoline 

P 

E 

F 

E 

F 

E 

T  etrahydrofur  an 

P 

E 

P 

E 

P 

E 

Toluene 

P 

E 

F 

E 

F 

E 

Xylenes 

P 

E 

F 

E 

F 

E 

Abbreviations  as  in  Table  SI. 

*  Note  that  for  analysis  NP-A2,  ratings  of  “Fail+”  were  treated  as  “Pass” 


Table  S3 

Additional  Solvents  and  Solvent  Mixtures  for  Test  PE-A2 


Solvent 

PE-A2 

P/F 

E/U 

Allyl  alcohol 

F 

E 

Nonylphenol 

P- 

E 

Acetonitrile  /  benzonitrile  87/13  v  /  v 

F 

E 

Acetonitrile  /  benzonitrile  77/23  v  /  v 

P- 

U 

Acetonitrile  /  benzonitrile  52/48  v  /  v 

P- 

E 

Acetonitrile  /  isopropanol  54  /  46  v  /  v 

F+ 

E 

Acetonitrile  /  n-butanol  53  /  47  v  /  v 

F 

E 

n-Butanol  /  propylamine  80  /  20  v  /  v 

F 

E 

Nitrobenzene  /  propylene  carbonate  52  /  48  v  /  v 

P- 

E 

Nitrobenzene  /  propylene  carbonate  22  /  78  v  /  v 

F 

E 

o-Dichlorobenzene  /  propylene  carbonate  60  /  40  v  /  v 

P 

E 

Quinoline  /  1,1,2,2,-Tetrabromoethane  52  /  48  v  /  v 

P 

E 

Abbreviations  as  in  Table  SI.  Note  that  the  expected  /  unexpected  behavior  for  the  other 
solvents  included  in  analysis  PE-A2  did  not  change  from  that  found  for  the  standard  analysis. 


The  following  Tables  (S4  -  S6)  provide  the  sources  and  values  for  molar  volume  and  cohesive 
energy  calculations  carried  out  using  the  “best  available  analog  group”  method.  The  method  is 
explained  in  detail  in  Appendices  A.l  (molar  volume)  and  A.2  (cohesive  energy)  in  the  main 
body  of  the  paper. 
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Table  S4 


Examples  Used  for  Estimating  the  Shrinkage  on  Dimerization 


Monomer 

Name 

Molar 
Volume 
(cc/mol) 
at  25  °C 

Source 

ID 

Dimer  Name 

Dimer 
Volume 
(cc/mol) 
at  25  °C 

Sourc 
e  ID 

Shrinkag 
e  on 

Dimer¬ 

ization 

(cc/mol)a 

Correctio 

n  to 

Monomer 

Volume 

(cc/mol)b 

n-Butane 

100.4 

lb 

n-Octane 

163.5 

lb 

37.3 

18.6 

n-Pentane 

116.2 

lb 

n-Decane 

195.9 

lb 

36.5 

18.2 

n-Hexane 

131.6 

la 

n-Dodecane 

228.6 

lb 

34.6 

17.3 

n-Heptane 

147.4 

lb 

n-Tetradecane 

261.3 

lb 

33.5 

16.8 

n-Octane 

163.5 

lb 

n-Hexadecane 

294.1 

lb 

32.9 

16.4 

Benzene 

89.4 

la 

Biphenyl 

155.1 

lb 

23.7 

11.8 

Isobutane 

98.1° 

2 

2, 5 -Dimethyl- 
hexane 

165 

3 

31 

16 

1-Chloro- 

butane 

105 

3 

1,8-Dichloro- 

octane 

178d 

4 

32 

16 

Notes:  a.  That  is,  twice  the  volume  of  the  monomer  minus  the  volume  of  the  dimer 


b.  That  is,  the  amount  subtracted  from  the  monomer  volume  used  to  estimate  the  volume 
of  the  attached  group,  e.g.  a  hexyl  group  would  be  17.3  cc/mol  smaller  than  hexane. 

c.  at  15  °C 

d.  at  20  °C 

Source  IDs: 

la.  Hansen,  C.  M.,  Hansen  Solubility  Parameters :  A  User's  Handbook.  2nd  ed.;  CRC  Press: 
Boca  Raton,  FL,  2007;  Table  A.l  (experimental  value) 

lb.  Hansen,  C.  M.,  Hansen  Solubility  Parameters:  A  User's  Handbook.  2nd  ed.;  CRC  Press: 
Boca  Raton,  FL,  2007;  Table  A.l  (calculated  value) 

2.  Air  Liquide  Product  Data  Sheet 

3.  Sigma- Aldrich  Product  Data  Sheet\ 

4.  Merck  Product  Data  Sheet 
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Table  S5 

Estimates  of  Molar  Volume  by  Best  Available  Analog  Group  Method 


Group 

Analog 

Analog 
Molar 
Volume 
(cc/mol) 
at  25  °C 

Source 

ID 

Correction 

Applied 

(cc/mol) 

Estimated 
Group 
Molar 
Volume 
(cc/mol)  at 

25  °C 

-Phenethyl 

Ethylbenzene 

123.1 

la 

17.5 

105.6 

-Styrenyl 

Styrene 

115.6 

la 

17.5 

98.1 

-Phenyl 

Benzene 

89.4 

la 

12 

77.4 

-1-Naphthyl 

Naphthalene 

111.5 

lb 

12 

99.5 

-Isobutyl 

Isobutane 

98. la 

2 

17.5 

80.6 

-Hexafluoro- 

isobutyl 

Hexafluoroisobutylene 

124 

3 

8.5b 

115.5 

-Trifluoro- 

propyl 

1,1,1  -Trifluoropropane 

96.9 

4 

17.5 

79.4 

Notes:  a:  at  15  °C 


b.  Based  on  applying  both  the  dimerization  correction  of  17.5  cc/mol,  and  an  additional 
correction  for  hydrogenation,  equal  to  the  difference  between  the  molar  volumes  of  isobutane 
(98.1  cc/mol  at  15  °C  from  source  2)  and  isobutylene  (89.4  cc/mol  at  25  °C  from  source  lb),  or  9 
cc/mol). 

Source  IDs: 

la.  Hansen,  C.  M.,  Hansen  Solubility  Parameters :  A  User's  Handbook.  2nd  ed.;  CRC  Press: 
Boca  Raton,  FL,  2007;  Table  A.l  (experimental  value) 

lb.  Hansen,  C.  M.,  Hansen  Solubility  Parameters:  A  User's  Handbook.  2nd  ed.;  CRC  Press: 
Boca  Raton,  FL,  2007;  Table  A.l  (calculated  value) 

2.  Air  Liquide  Product  Data  Sheet 

3.  DuPont  Product  Data  Sheet 

4.  Mattox,  D.  M.  The  Foundations  of  Vacuum  Coating  Technology ;  Noyes  /  William  Andrew 
Publishing:  Norwich,  NY,  2003. 
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Table  S6 

Estimates  of  Cohesive  Energy  by  Best  Available  Analog  Group  Method 


Group 

Analog 

Analog  Cohesive  Energy 
(J/mol)  at  25  °C 

Source 

ID 

Dispersive 

Polar 

H-Bonding 

-Phenethyl 

Ethylbenzene 

39000 

40 

240 

la 

-Styrenyl 

Styrene 

40000 

120 

1900 

la 

-Phenyl 

Benzene 

30000 

0 

360 

la 

-1-Naphthyl 

Naphthalene 

41000 

450 

3900 

lb 

-Isobutyl 

Isobutane 

19000a 

0 

0 

2 

-Hexafluoroisobutyl 

Hexafluoroisobutylene 

22000b 

0 

0 

3 

-T  rifluoropropyl 

1,1,1  -Trifluoropropane 

13000° 

3900d 

0d 

lb, 4, 5 

Notes:  a.  based  on  AHvap  at  25  °C 

b.  Based  on  AHvap  at  25  °C  determined  from  the  Clausius-Clapeyron  fonn  of  the  listed  vapor 
pressure  data,  assumed  entirely  dispersive,  and  a  correction  for  hydrogenation,  equal  to  the 
difference  between  the  dispersive  cohesive  energies  of  isobutane  (1,9200  J/mol  at  15  °C  from 
source  2)  and  isobutylene  (18,800  J/mol  at  25  °C  from  source  lb),  or  400  J/mol). 

c.  Based  on  the  total  (from  AHvap  at  25  °C  reported  in  source  4)  minus  the  expected  polar  and 
hydrogen-bonding  contributions 

d.  The  values  for  the  polar  and  H-bonding  cohesive  energies  are  taken  from  those  reported  in 
source  lb  for  1,1,1-trifluoroethane,  corrected  for  temperature  by  scaling  the  computed  values  of 
EC0h  (based  on  AHvap)  at  an  assumed  174  K  (from  source  lb)  to  that  at  298  K  (from  source  5), 
Note  that  the  values  reported  in  source  5  at  174  K  match  those  reported  in  source  lb,  both  for 
molar  volume  and  AHvap.  These  values  are  assumed  not  to  be  altered  by  addition  of  a  methylene 
group. 

Further  explanation: 

According  to  source  lb,  Ecoh  for  1,1,1-trifluoroethane  is  21,200  J/mol,  with  ED  =  13,800 
J/mol  and  EP  =  7,400  J/mol.  Accounting  for  PV  yields  AHvap  =  22,600  J/mol.  According  to 
source  5,  AHvap  at  174  K  equals  approximately  this  value,  however,  at  298  K  (making  use  of  the 
thermodynamic  tables),  AHvap  =  13,700  J/mol,  or,  equivalently,  Ecoh  =  11,300  J/mol.  Hence,  the 
component  values  from  source  lb  have  been  scaled  by  the  factor  (1 1,300)  /  (21,200),  resulting  in 
estimates  of  Ed  =  7,400  J/mol  and  Ep  =  3,900  J/mol  at  298  K  for  1,1,1-trifluoroethane. 

For  1,1,1-trifluoropropane,  the  value  of  AHvap  at  298  K  from  source  4  is  19,600  J/mol,  or 
Ecoh  =  17,100  J/mol  after  accounting  for  PV.  In  order  to  estimate  components,  it  is  assumed  that 
the  dipoles  in  1,1,1-trifluoropropane  are  simply  diluted  versions  of  the  dipoles  in  1,1,1- 
trifluoroethane,  hence  EP  is  assumed  equal  for  both  compounds,  at  3,900  J/mol.  It  is  assumed 
that  Eh  also  remains  zero,  hence  the  remaining  13,200  J/mol  is  assigned  to  ED. 

Source  IDs: 

la.  Hansen,  C.  M.,  Hansen  Solubility  Parameters :  A  User's  Handbook.  2nd  ed.;  CRC  Press: 
Boca  Raton,  FL,  2007;  Table  A.l  (experimental  value) 

lb.  Hansen,  C.  M.,  Hansen  Solubility  Parameters:  A  User's  Handbook.  2nd  ed.;  CRC  Press: 
Boca  Raton,  FL,  2007;  Table  A.l  (calculated  value) 
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2.  CRC  Handbook  of  Chemistry  and  Physics,  82nd  ed.  (D.  R.  Lide,  ed.);  CRC  Press:  Boca 
Raton,  FL,  2001;  p.  6-1 13. 

3.  DuPont  Product  Data  Sheet 

4.  Mattox,  D.  M.  The  Foundations  of  Vacuum  Coating  Technology ;  Noyes  /  William  Andrew 
Publishing:  Norwich,  NY,  2003. 

5.  NIST/TRC  Web  Thermo  Tables  (WTT);  available  online  at 

http://webbook.nist.gov/cgi/fluid.cgi?ID=C420462&Action=Page,  with  additional  information 
available  at  http://wtt-lite.nist. gov/wtt-lite/index.html?cmp=l.l.l-trifluoroethane,  accessed 
March  1,  2012. 
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